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ABSTRACT 

Because  materials  are  often  stressed  into  the  clastic  ranpe  in 
light  weight  construction,  a need  exists  for  information  concerning 
the  plastic  behavior  of  structural  members.  A plate  subjected  to  lateral 
loads  is  such  a structural  member.  In  this  report,  we  use  several 
general  methods  of  plasticity  to  obtain  information  about  the  plastic 
behavior  of  circular  plates  and  we  suggest  how  these  methods  cun  be 
extended  to  rectangular  plates. 

In  this  report,  the  kinematic  relations  connecting  the  displacements 
of  the  neutral  surface  to  the  strains  are  those  of  the  small  deflection 
theory  of  bending.  The  stress-strain  relations  are  those  of  the  theory 
of  plastic  deformations  (secant  modulus  theory).  In  the  last  chapter, 
the  comparison  of  the  theory  of  plastic  flow  and  the  theory  of  plastic 
deformations  indicates  that  only  slight  differences  would  be  fourid  in 
numerical  results  calculated  by  either  of  these  two  theories. 

The  four  approximate  methods  which  are  applied  to  circular  plates 

are 


1) .  Sokolovsky's  Method, 

2) .  Iteration  Method, 

3) .  Potential  Energy  Method,- 

U).  Complementary  Potential  Energy  Method. 

Sokolovsky  reduces  the  equations  relating  the  moments,  curvatures,  and 
loads  of  the  plate  to  two  simultaneous  first  order  non-linear  differential 
equations.  These  equations  are  then  solved  by  numerical  integration.  The 
iteration  method,  developed  by  Ilyushin,  is  adapted  for  the  bending  of 
circular  plates.  In  this  method  we  separate  a non-linear  second  order 
differential  equation  into  a linear  portion  and  a non-linear  portion.  The 
effect  of  the  non-linear  portion  on  the  solution  is  found  by  the  iteration 
procedure.  The  principle  of  minimum  potential  energy  is  used  to  estimate 
the  circumferential  cui vature  of  the  plate  by  means  of  functions  containing 
arbitrary  parameters.  These  parameters  are  evaluated  by  Galerkin's  Method. 
The  principle  of  minimum  complementary  potential  energy  is  used  to  estimate 
the  radial  bending  moment  in  the  plate  by  a similar  process. 

Numerical  results  are  obtained  by  all  four  methods  for  a simply 
supported  uniformly  loaded  plate  for  a material  where  the  second  stress 
invariant  S i3  related  to  the  second  strain  invariant  E by 

S - 2G(l  - XE2)E. 

A comparison  of  the  results  of  the  above  methods  is  riven  in  the  last 
chapter  of  this  report.  Also,  in  this  chapter,  the  last  three  of  the 
above  approximate  methods  are  developed  for  rectangular  plates,  but  no 
numerical  calculations  are  made. 
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1.1  The  Problem  Before  Us. 

^he  circular  plate  symmetrically  loaded  by  lateral  forces  is  a 
structural  member  which  the  practicing  engineer  often  encounters.  A 
schematic  diagram  of  the  problem  is  shown  in  figure  1.  The  circular 
symmetry  of  the  loading  enables  the  deflection  and  stress  condition  of 
the  plate  to  be  represented  as  functions  dependent  upon  the  radial 
distance  r.  Thus,  the  physical  relationships  are  expressed  in  terras 
of  ordinary  differential  equations. 

This  work  considers  the  small  deflection  theory  of  plates  due 
to  lateral  loads.  The  strains  of  the  neutral  surface  and  the  quadratic 
terras  in  the  strains  due  to  bending  are  neglected.  The  laws  of  the 
theory  of  plastic  deformations  are  used  instead  of  elastic  stress- 
strain  relationships. 

Naghdi  (l)  considers  the  large  deflection  theory  of  circular  plates 
with  stretching  of  the  neutral  surface.  Oleyxal  (2),  Mostov  (3),  and 
Hill  (U),  approximate  the  stress  condition  for  the  large  deflections 
of  thin  plates  by  considering  only  membrane  stresses  and  strains. 

1.2  Objective 

For  a general  problem,  such  as  the  bending  of  a plate,  several 
procedures  are  available  for  obtaining  a solution.  The  preferable 
method,  in  any  particular  case,  depends  upon  the  accuracy  and  infor- 
mation desired.  The  mo3t  desirable  method  is  one  which  yields  the 
exact  solution.  For  the  plastic  bending  of  a circular  plate  no  such 
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exact  solution  has  been  found.  Therefore,  this  thesis  applied  several 
general  methods  of  plasticity  to  obtain  approximate  solutions  and  a 


comparison  is  made  of  the  accuracy  and  labor  involved  in  these  methods. 

In  the  actual  computation  of  the  example  problem,  a dimensionless 
parameter  appears  in  terras  of  which  the  range  and  accuracy  of  the 
various  methods  are  indicated.  The  comparison  of  the  results,  obtained 
by  the  different  methods  of  solution  applied  to  the  problem  of  the 
circular  plate,  should  give  some  insight  into  those  problems  which 
present  greater  mathematical  difficulties,  such  as  the  bending  of 
rectangular  plates. 

1.3  The  Kinematic  and  Equilibrium  Equations  for  the  Symmetrical  Bend- 
ing of  Circular  Plates. 

This  section  states  the  kinematic  and  equilibrium  conditions  for 
the  symmetrical  bending  of  circular  plates.  A complete  presentation 
of  what  is  mentioned  here  is  given  by  Timoshenko  (5>)  (Chap.  3). 

The  following  assumptions  are  made: 

a) .  The  defloction  of  the  neutral  surface  of  the  plate  is  of 
such  a magnitude  that  the  quadratic  terras  due'to  bending  are  neglected 
in  the  strain  formulae. 

b) .  The  elongation  of  the  neutral  surface  and  the  membrane  stresses 
are  neglected. 

c) .  Normals  to  the  neutral  surface  before  deflection  are  normal 
to  the  qputral  surface  in  the  deflected  state. 

d) .  The  normal  3tress  tr  ie  assumed  zero  throughout  the  plate. 

z 

e) .  One  additional  condition  is  assumed  due  to  the  nature  of 

the  plastic  laws.  This  condition  is  that  the  loading  is  applied  proportionally 

i 
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and  increases  monotonously.  Hence,  only  loading  is  considered. 

In  stating  the  kinematic  relations  of  the  plate  the  following 
notations  are  used. 

w - deflection  of  the  neutral  surface  from  the  unstrained  to  the 
strained  condition. 

Oq  ■ l/Pg  “ circumferential  curvature. 

a • l/p  * radial  curvature 
r /J  r 

u “ radial  displacement 
z ■ distance  from  the  neutral  surface. 

The  curvatures  are  expressed  in  terms  of  the  deflection  as 


/ dw 

a9 ' - a? 


and  a 


(1.1) 


The  geometry  of  figure  2 suggests  the  physical  significance  of  aQ 
and  a . The  radial  radius  of  curvature,  -P , is  the  radius  of  curvature 
of  the  intersection  of  the  neutral  surface  and  the  vertical  plane  con- 
taining a radius  of  the  plate.  The  circumferential  radius  of  curvature, 
-P  , is  the  radius  of  curvature  of  the  intersection  of  the  neutral 

v 

surface  and  the  plane  perpendicular  to  the  tanrent  of  the  curve  in  the 
neutral  surface  at  9 ■ const. 

The  geometry  of  figure  3 gives 

dw 

U " " Z 3r  * 

Due  to  the  circular  symmetry  of  the  problem,  the  circumferential  dis- 
placement is  zero. 

From  known  formulae  for  the  strains  in  terms  of  displacements,  the 

equations  for  the  radial  strain,  e , and  the  circumferential  strain,  e , 

r y 


i 


u 


© r 


,d2w 

d7 


z dw 

r 3r  M© 


(1.2) 


Another  important  equality  derived  from  (l.l)  i3 
d(ra©) 


(1.3) 


The  following  notation  is  used  to  fonailate  the  equilibrium 
equations  x 

p(r)  ■ lateral  load  per  unit  area. 

Q ■ radial  shearing  force  per  unit  length, 
r 

M ■ radial  moment  per  unit  length, 
r 

Mg  ■ circumferential  moment  per  unit  length. 

A pictorial  representation  of  these  quantities  is  given  in  figure  U. 
The  definition  of  the  moments  is 


Mr  " 2 \ °~r  Z dz’ 


(1.U) 


*6  " 2 } 


where  0~  is  the  radial  stress  and  <Ti  is  the  circumferential  stress, 
r 9 

It  has  been  assumed  that  the  radial  and  circumferential  stresses  have 
point  symmetry  about  z ■ 0. 

Consideration  of  the  equilibrium  of  the  forces  on  the  circular 
section  of  the  plate  shown  in  figure  Ujl  yields 


rQr  ■ J prdr. 


(1.5) 


i 
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The  equilibrium  equation  is  found  by  using  the  above  equality  and  the 
equilibrium  of  the  moments  on  the  element  shown  in  figure  Ub 


d(rMr) 

3r 


5% 


ordr  ■ 0. 


(1.6) 


1.U  Development  of  the  Stress-Strain  RelationshiDS. 

This  work  is  not  intended  to  be  a discussion  of  the  relative 

merits  of  the  different  theories  of  plasticity;  therefore,  this  section 

briefly  presents  the  laws  of  plasticity  used  in  this  paper.  A complete 

development  of  the  relationships  used  in  this  section  is  given  by 

Sokolovsky  (6)  (chap.  1).  A thorough  explanation  of  the  yield  condition 

is  found  in  the  book  of  Hill  (7)  (chap.  2).  The  theory  of  plastic 

deformations  (secant  modulus  theory)  is  used. 

The  plastic  stress-strain  relationships  are  now  introduced.  The 

principal  stresses  are  <y-  , cra>  and  o~  . The  principal  strains  are 

r y z 

e , e_,  and  e . 
r*  9’  z 


The  mean  stress  is 

1/3  (cr  + a; g ♦ c r ). 

The  mean  strain  is 

e - 1/3  (er  ♦ eQ  ♦ ez). 

The  stress  deviations  are 
s ^ ■ cr^  -<r(i  “ r,  9,  z). 

The  strain  deviations  are 
®i  ■ «i  - * (i  - r,  0,  z). 

The  second  stress  invariant  is 


(1.7) 


(1.8) 


(1.9) 


(1.10) 


S'  “ I (3r2  + s02  * O' 


(1.11) 
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The  second  strain  invariant  is 


E'  “ I (6r2  * e02  * ez2>‘ 


(1.12) 


The  following  assumptions  of  the  secant  odulus  theory  are  used, 

a).  The  stress  deviations  are  nroportional  to  the  correspond inf- 
strain  deviations. 


S 3.  3 

r _ 6 _ z 

e e0  e 

r 0 z 


With  the  use  of  alpebraic  operations  and  the  definitions  (1.11)  and 
(1.12),  the  above  relation  i3  formulated  as 


S 3 _ 3 q 

e ert  e B 

r 0 z 


(1.13) 


b) .  The  second  stress  invariant  is  a function  of  the  second  strain 
invariant.  This  relationship  is  written  as 

8-20[l-  f (E)j  E.  (l.Ui) 

If  the  relationship  between  S and  E has  point  symmetry  about  the  origin, 
thenf(E)  is  an  even  function  of  E.  Thus,  the  function,  f(E),  can  be 
expressed  as 

f (E)  - £ XnE2n,  (1.15) 

s=f  n 

where  the  X's  are  constants. 

c) .  The  material  is  incompressible, 

• ■ J (er  * ‘e  * •,>  ■ 0 
from  which  it  follows  that 

*,  ■ - (e6  * er)* 


(1.16) 


This  assumption  greatly  simplifies  the  .lathematical  formulation  of  the 
problem. 

The  principal  stress  deviations  and  stress  invariant  are  obtained 

by  using  equations  (1.7),  (1.9),  (1.11),  and  assumption  d in  section  1.3 

that  cr  “0* 
z 


8r  " f (<rr  " ~T')’  S6  " | (cr0  " -£)>  \ m - \ {crr  * °9)  (1*17) 


S2  - ^ (c rr2  - arc rQ  ♦ <rQ2).  (1.18) 

The  principal  strain  deviations  and  strain  invariant  are  obtained  with 
the  aid  of  equations  (1.8),  (1.10),  (1.12),  and  (1.16): 

er  - V e8  " *9'  - - ' (er  * *9>  (1-19) 

2 


E2-e2+ee_+e_ 
r r 0 0 


(1.20) 


With  the  expressions  (1.13),  (1.15),  (1.17),  (1.19),  and  (1.20),  the 

principal  stresses  cr  and  cr  are  expressed  in  terms  of  the  principal 

I*  w 


strains : 


or  - UG 


°e  • W 


t v 1 ■»  / 2 2 <n  / 0 

1 - 1 , X (e  ♦ ® ) («  ♦ -x— 

n r r 0 0 7 ' r 2 

n«l 


1 -gXn  (er2  * ere9  * e 9 


"j  (te  * T-' 


(1.21) 


These  are  the  stress-strain  relations. 

1.5  Formulation  of  the  Basic  Plate  Equations. 

This  section  Dresents  the  formulae  which  are  used  in  the  solution 


of  the  bending  of  a circular  plate 


8 


By  substituting  exoressiens  (1.1)  and  (1.21)  into  equations  (l.u) 
and  integrating,  we  find  the  bending  noments 


Mr  - D (1  - B)  (or  - -|-) 


^ - D (1  - B)  (o9  ♦ -|-), 


where 


D - G 


(1.22) 


(1.23) 


The  integration  of  (l.U)  is  possible  because  the  a and  aQ  are  constant 

r w 

for  a riven  r by  assumption  c in  section  1.3.  The  constant  D is  the 
customary  expression  for  the  flexural  rigidity  of  a plate  in  elasticity. 
The  quantity  B expresses  the  non-linearity  characteristic  of  the  plastic 
stress-strain  relationship. 

The  solution  for  the  bending  of  the  ciicula.'  plate  is  found  when 
M^,  Mg,  «r,  Og,  and  w are  found  as  functions  of  r.  The  five  relation- 
ships in  (1.1),  (1.3),  (1.6),  (1.22)  are  sufficient  for  the  determi- 
nation of  the  above  functions.  The  combination  of  these  five  equations 
yields 


i 


d3w  1 d2w  1 dw  1 fr 

7T  +F72  -TdF"TJF  ) prdr 

dr  dr  r o 


(1.2U) 


d ff^L  . ^ 3Xnh  n /d2w)2  /dfwVl  dwV  /l  dw)2  n( 

* dr((dr2  2r  £=1  2n+3  l d 7/  tdr2/Vr  ^ V**  drJ  j 

1 1 d2v  1 dv")  3Xn^  * ( d2w)2  /d2w)/l  dw)  ^/l  dw)2  n 

2l?^‘7“J?JS  WA72?A?a?JJ  . 


The  values  M , MQ,  a , aQ,  and  w are  not  found  by  solvinr  equation 
r w r y 

(l.2ii)  directly.  Different  dependent  variables  are  introduced  in  this 
thesis  which  reduce  the  numerical  work  in  the  various  approximate 
methods. 

Table  I is  compiled  for  future  reference.  The  formulae  in  this 
table  are  (1.1),  (1.3),  (1.6),  (1.22),  and  (1.23). 


0 r dr 

d(ro  ) 
x - — 

r dr 


TABLE  I 


(II.) 


d(rMr) 


-vf 


prdr  B 0 


a a 

Mr  - D(1  - B)  (ar  ♦ -J-),  M9  - D(1  - B)  (ae  * ^-) 


(III.) 


(IV.) 


OO  3A  h 


-jfey  («/  . araQ  . aQAr 


! 


CHAPTHl  II 


THE  METHOD  OF  SOKOLOVSKY 

2.1  Mathematical  Formulation  of  the  Problem. 

Sokolovsky  (6)  (chap,  ill)  treats  the  problem  of  the  bending  of 
a circular  plate  symmetrically  loaded.  In  orinciDle,  Sokolovsky 
reduces  formulae  II,  III,  and  IV  into  two  simultaneous  first  order 
ordinary  differential  equations.  By  the  introduction  of  two  new 
variables,  a magnitude  factor  A and  an  angular  factor  to,  which  re- 
place the  curvatures  and  o0,  these  two  differential  equations  can 
be  expressed  in  a convenient  form.  The  variables  A and  co  are  defined 

in  terms  of  a and  as 
r 6 

sin  (u  - ti/6),  aQ  * sin  (a>  + n/6).  (2.1) 

With  the  use  of  (2.1),  the  following  identity  for  A i3  formed: 

A2  - a/  ♦ ora0  + aQ2.  (2.2) 

A combination  of  (1.2),  (1.20),  and  (2.2)  yields 

B2  - z2A2,  (2.3) 

which  shows  that  the  factor  z is  the  proportionality  constant  between 
E and  A.  This  same  proportionality  exists  between  the  principal  strains 
and  the  curvatures  as  seen  in  (1.2). 

A graphical  representation  of  a>  is  given  in  figures  5a  and  5b. 

The  normal  of  the  plane  of  figure  5b  makes  equal  angles  with  the 
three  orthogonal  axes  ar’,  a0',  - (eM  ♦ a0')  of  figure  5a.  This 
plane  is  called  the  octahedral  plane.  The  axes  shown  in  figure  5b 
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11 


are  the  projections  of  the  axes  of  figure  5a.  A vector  with  the 
components  \/2/3  ar' > /2/3  a$' > /2/3  |j-  (ttr'  + aQ,)~ |>  has  bhe 

length  2//3A  and  lies  in  the  plane  of  figure  5b.  The  vector  2//3A 
lies  at  an  angle  co,  as  shown  in  figure  5b.  It  has  the  orthogonal 
projections  a r and  on  the  axes  a and  in  the  octahedral  plane. 
When  the  proportional  nrincipal  strains  aie  substituted  for  the 
curvatures  in  the  a’^ove  graphical  representation,  the  vector  2//3A 
becomes  the  octahedral  shearing  strain. 

With  the  substitution  of  (2.1)  and  (2.2)  in  IV  and  V of  Table  I, 
the  moments  are  expressed  as 

■ D (1  - B)  A cos  (co  + n/6),  (2.1i) 

Mg  - D (1  - B)  A cos  (co  - n/6). 


where 


B - 


3Xn 

2n  + 1 


(hA) 


2n 


(2.5) 


The  substitution  of  (2.1)  into  II  of  Table  I yields 

r jjsin  (co  + n/6)  gp  + A cos  (co  + n/6)  g^j  ♦ /3A  sin  co  ■ 0.  (2.6) 
Similarly,  the  substitution  of  (2.h)  into  II  of  Table  I gives 
r jj:os  (co  n/6)  (1-  gy)  gp  - (1  - B)  A sin  (co  + n/6)  g^j 

-f 


- (1  - B)  A sin 


rp  (r)  dr/D  ■ 0. 


(2.7) 


These  two  equations  are  reorganized  into  a more  convenient  form  for 
numerical  calculations  when  the  derivatives  of  co  and  A are  separated. 


This  separation  yields 
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» 


(l-B)-A-jy  + A-jy  sin^  (co  ♦ n/6) 


gp  - A^  3 in  u cos  (co  ♦ n/6) 


♦ A (1  - B)  sin  2 £o  - sin  (co  ♦ n/6)  ^ rp(r)  cir/D  ■ 0 


(2.8) 


(1  - 


. dB  . dB  , 2 , 
B)  " A dA  + A dA  Sin  ( 


“ ♦ n/6)  -j| 


♦ 2A(1  - B)  sin 


CO 


♦ cos  (to  ♦ n/6)  ^ rp(r)  dr/D  - 0. 


(2.9) 


Equations  (2.8)  and  (2.9)  are  the  desired  first  order  ordinary 
differential  equations.  They  can  be  solved  by  numerical  methods  when 
definite  stress-strain  laws  and  boundary  conditions  are  given.  When 
r ■ 0 and  A / 0,  equations  (2.8)  and  (2.9)  both  yield  co  * 0 and  A 
is  finite.  The  substitution  of  co  * 0 into  (2.1*)  gives 

Mr  “ % " ^ D(1  " B)  A* 

Expression  III  Table  I also  yields  that  Mr  * M at  the  center  of  the 
plate  where  r * 0.  Thus,  one  boundary  condition  is: 

At  the  center  of  the  plate  where  r - 0 

M - M or  co  - 0 . (2.10) 

r 0 

Two  boundary  conditions  are  necessary  to  solve  equations  (2.8)  and 
(2.9) • The  second  boundary  condition  is  obtained  by  physical  restrictions 
on  tho  curvatures  or  moments  at  the  outer  edge  of  the  plate.  With  the 
use  of  either  equations  (2.1)  or  (2.1*)  and  the  given  boundary  conditions 
at  the  outer  edge  of  the  plate,  values  for  A or  co  are  determined  for 
this  outer  boundary.  In  the  next  two  sections,  a method  of  solution 
is  demonstrated  by  solving  an  example  problem. 
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2.2  The  Illustrative  Example. 

An  illustrative  example  of  a simply  supported  plate  with  constant 
pressure  p is  solved  to  indicate  the  method  of  solution  for  (2.8)  and 
(2.9).  The  stress-strain  relation  (l.lU)  is  assumed  as 

S - 2G  (1  - XE2)  E . (2.11) 

A graphical  representation  of  this  curve  is  given  in  figure  6 which 
shows  a horizontal  tangent  at  E ■ Eq.  From  the  condition  of  a hori- 
zontal tangent  at  Eq,  the  constant  X is  determined  by 

d S/dE  - 2G  (1  - 3XEq2)  - 0. 


from  which  it  follows 


X - 1/3  Eq  . 


(2.12) 


Two  new  dimensionless  variables  a and  y are  introduced  to  simplify 
the  form  of  (2.8)  and  (2.9)  for  this  example.  These  variables  are 
expressed  in  terms  of  A and  r as 

a - A/Aq,  y - (Xl/2  h p/D)l/2r,  (2.13) 


where  the  constant  A is  defined  by 

o J 


h2A  2 - 1/3  X. 


(2.Ui) 


The  maximum  E must  occur  at  z ■ h,  and  Eq  is  the  maxinum  permissible 

strain  invariant.  A combination  of  (2.3),  (2.12),  and  (2.lU)  gives 
2 2 2 

Eq  • h Aq  j thus,  Aq  can  be  interpreted  as  a measure  of  the  maximum 
permissible  strain  condition  in  the  plate.  This  interpretation  limits 
the  value  of  a to  0— a £1.  The  factor  B,  which  is  given  in  expression 
(2.5),  is  given  in  terms  of  a as 


Ill 


B - | Xh2A2  -0.2a2  (2.15) 

for  the  stress-strain  relation  (2.11).  When  the  pressure  p is  constant, 
the  load  integral  is 


(2.16) 


o 

With  the  use  of  (2.15)  and  (2.16),  equations  (2.8)  and  (2.9)  are 
expressed  in  terms  of  the  variables  co,  a,  and  y as 

ya  Jl  - 0.6a2  ♦ O.Ua2  sin2  (co  ♦ n/6)J  — - 0.6928  a^  sin  co  cos  (co  ♦ n/6) 
♦ a (1  - 0.2a2)  sin  2 co  - 0.866  y2  sin  (co  ♦ n/6)  ■ 0 (2.17) 


y jj.  - 0.6a2  ♦ O.Ua2  sin2(co  + «/6~)~|  ^ ♦ 2a  (1  - 0.2a2)  sin2co 

♦ 0.866  y2  cos  (co  ♦ n/6)  * 0 (2.18) 

For  a simply  supported  plate  the  radial  moment  is  zero  at  the 
outer  edge  of  the  plate.  From  (2.U),  this  boundary  condition  is 
evaluated  as : 

1) .  At  the  outer  edge  of  the  plate,  r * b 

Mr  - 0 or  co(T)  - n/3  where  I - [x1/2  h pb2/oj  , (2.19) 

2) .  At  the  center  of  the  plate,  r ■ 0,  expression  (2.10)  gives 

- Mq  or  co  - 0 at  y - 0.  (2.20) 

Conditions  (2.19)  and  (2.20)  determine  unique  solutions  of  (2.17) 
and  (2.18). 

A step  by  step  numerical  integration  of  (2.17)  and  (2.18)  can  be 
started  when  a value  for  *a"  at  y ■ 0 is  chosen.  This  step  by  step 
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I . 

I 


integration  is  continued  until  u * n/3,  and  by  condition  (2.19)  the 
value  of  Y,  where  u>  ■ n/3,  determines  the  outer  boundary  of  the  plate. 
When  Y is  orescribed,  a series  of  starting  values  for  •a"  are  assumed. 
These  starting  values  are  adjusted  until  one  fives  the  correct  Y at 
W(Y)  ””  1°  the  next  section  the  numerical  procedure  is  explained 

by  an  example  where  the  initial  conditions  are  y * 0,  u>  m 0,  and 
a - 0.5. 

2.3  The  Numerical  Integration: 

A numerical  method  used  for  computational  work  should  have  a 
checking  procedure.  When  evaluation  of  the  numerical  steps  )>econes 
tedious,  calculation  errors  arise  unless  an  automatic  check  is  available. 

An  accurate  procedure  satisfying  the  above  conditions  is  Milne's 
Method.  Scarborough  (8)  (pp  2lt5  and  295)  gives  a thorough  explanation 
of  this  process.  In  thi§  method,  the  following  formulae  are  used: 

Milne's  formula 

^1  " “n-3  + (2  3?|n-2  " 3?  In-1  * 2 3? In)  » (2*2l) 

Simpson's  formula,  which  is  used  as  a check, 

“nil  " Vl  + H (af |n-l  + ^afln  * afln+l)  (2‘22) 

and  the  error  term 

error  - ^ - J^)  . (2.23) 

Milne's  fornula  passes  a curve  through  four  points  and  estimates 
the  value  of  the  function  at  the  fifth  point  by  an  extrapolation  of 
the  curve,  (see  figure  7).  When  the  value  of  is  calculated  with 
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the  aid  of  formula  (2.21),  the  given  differential  equation  is  used 


dco 


to  find  the  slope  |n+l*  Ane  va-LUC  D1  | n+1 


The  value  of 


dco 


obtained  from  the 


above  procedure  is  used  in  Simpson’s  formula,  (2.22),  to  make  a 

(2) 

second  approximation  for  co^+'. 

When  exact  values  are  used  in  (2.21)  and  (2.22),  formula  (2.21) 
gives  an  which  is  too  small  while  (2.22)  gives  an  which 
is  too  large.  A consideration  of  the  terras  which  are  neglected  in 
the  extrapolation  formilae  (2.21)  and  (2.22)  shows  that  the  maximum 
deviation  from  the  true  curve  is  given  by  the  formula  (2.23).  Thus, 
formula  (2.22)  provides  a check  on  the  numerical  calculations  for 
(2.21),  and  (2.23)  gives  a means  of  determining  the  accuracy  of  the 
approximation  method. 

Before  Milne's  Method  can  be  applied,  the  values  of  the  function 
and  its  first  derivatives  at  four  points  must  be  obtained.  The  value 
of  the  function  at  the  second  point  is  found  with  the  aid  of  Euler's 
averaging  method.  The  formula  used  in  this  method  is 


co  ♦ 
o 


dcol 


, dco 

o 3y 


0 * 


(2.2U) 


The  value  of  the  function  at  the  third  and  fourth  points  is  found  by 
applyinp  Simpson's  formula,  (2.22).  The  numerical  procedure  is 
explained  by  showing  the  steps  involved  in  solving  equations  (2.17) 
and  (2.18)  when  the  initial  conditions  are  co^  ■ 0 and  a^  * j. 

The  solutions  of  equations  (2.15)  and  (2.16)  for  dco/dy  and 
da/dy  are 


— r-1 
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dto/dy 


- da/dy 


^.69?B  a^sin  to  co3  (to  ♦ n/6 ) - a(l  - 0.2a2 )sin  2 to 

♦ 0.866  y2sin  (to  + jtfl 
ya 


JjL  - 0.6a2  ♦ O.lia2  sin"’  (to  ♦ n/6}] 


2a(l  - 0.2a2)  sin2to  ♦ 0.866  y2  cos  (to  ♦ n/6) 


(2.25) 


y^  - 0.6a2  ♦ o.ija2  si n2  (to  ♦ n/6)j  (2.26) 


For  the  initial  starting  value  of  to  ■ 0,  formulae  (2.25)  and  (2.26) 

yield  indeterminate  equations  for  g—  and  These  equations  are  of 

the  form  0/0.  The  evaluation  of  (2.25)  and  (2.26)  at  to  ■ 0 can  be 

2 

accomplished  by  expanding  sin  to  and  sin  to  about  to  " 0 in  terms  of  y. 
Theae  expansions  are 


, dxo 

sin  to  - y ^ 


sin  2 to  - 2y  ^ 


d2to| 

77  < 

dy  |o 

2 d2to 
dy 


With  the  substitution  of  these  expansions  into  (2.25)  and  (2.26), 
the  limits  of  the  expressions  for  the  derivatives  as  y-»0  are 


lira 


dto  ^ . dto 

y-»  0 3y  3y 


lim  da  m 


|o  * y-*0  3y  ® ’ 

This  procedure  yields  starting  values  for  dto/dy  and  da/dy  of 
dto 


3y  ’ ° » 3y  “ ° * (2‘27) 

For  an  increment  in  y,  steps  of  0.05  are  chosen.  The  initial  increment 
in  y should  be  small  because  of  the  irregular  behavior  of  to  at  the 
origin.  Also,  the  checking  methods  do  not  apply  until  the  values  for 
a,  to,  dto/dy,  and  da/dy  are  calculated  for  four  values  of  y. 


i 


18 


First  step:  Calculation  for  the  point  y ■ 0.0?. 

The  values  co  ■ 0 and  a ■ 0.5  in  the  row  y * 0 in  Table  II  are 
the  prescribed  initial  conditions.  First  approximations  for  co  and  a at 
y ■ 0.0$  are  assumed  to  be  co  ■ 0 and  a * 0.5.  First  aoproximations 
for  dco/dy  and  da/dy  at  y * 0.05  are  obtained  by  using  these  values  of 
co  and  a in  formulae  (2.25)  and  (2.26).  The  newly  calculated  dco/dy 
and  da/dy  are  substituted  in  formulae  (2.2U)  to  obtain  second  approxi- 
mations of  co  and  a.  These  values  of  co  and  a can  be  substituted  in 
formulae  (2.25)  and  (2.26)  to  calculate  second  approximations  for 
dco/dy  and  da/dy.  With  these  values  of  the  derivatives,  third  approxi- 
mations for  co  and  a are  obtained  from  the  expression  (2.2U). 

This  process  of  approximation  is  repeated  until  the  values  for 
dco/dy  and  da/dy  remain  stationary  in  which  case  the  co  and  a given  by 
(2.2U)  remain  constant.  These  stationary  values  of  co,  "a"  and  their 
derivatives  are  given  in  the  row  y * 0.05  ' Table  II. 

In  the  above  procedure,  the  values  of  co  oscillate  about  the 
stationary  co  with  a decreasing  amplitude.  With  the  use  of  slight 
alterations  in  the  value  of  co  towards  the  center  of  the  oscillation, 
the  number  of  calculations  can  be  reduced.  The  magnitude  of  these 
alterations  can  be  estimated  after  the  first  two  calculations. 

Second  Step:  Calculations  for  the  points  y “0.10  and  y * 0.15. 

Approximate  values  of  co  and  a at  y ■ 0.10  are  obtained  with  the 
use  of  Simpson’s  formula,  (2.27).  In  Simpson’s  formula,  however,  the 
value  of  the  derivative  at  the  point  to  be  calculated  must  be  known. 
Therefore,  an  estimation  of  this  derivative  is  made  by  a linear  extrapolation 


i 
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of  the  two  preceding  derivatives.  The  approximate  values  of  u and  a 
calculated  by  Simpson’s  formula  are  substituted  into  expressions  (2.25) 
and  (2.26)  to  obtain  new  estimates  of  the  derivatives.  These  values 
of  the  derivatives  are  used  to  make  new  approximations  for  w and  "a" 
with  Simpson’s  formula.  Two  calculations  are  usually  sufficient  to 
obtain  the  desired  accuracy  for  the  value  of  the  functions  and  their 
derivatives.  Exactly  this  same  procedure  is  used  to  make  the  calculations 
for  a and  "a"  at  y ■ 0.15.  The  results  of  these  calculations  are  shown 
in  Table  II. 

Third  Step:  Calculations  for  further  values  of  y. 

Once  four  values  of  <0,  a,  and  their  derivatives  are  obtained,  the 
regular  numerical  integration  procedure  is  adopted.  Formula  (2.21) 
gives  the  values  and  a^|  for  y • 0.2.  These  values  are  substi- 

tuted into  formulae  (2.25)  and  (2.26)  to  calculate  da>/dy  and 
da/dy  These  values  of  co,  a,  and  their  derivatives  are  entered 

in  the  row  y - 0.2. 

All  the  quantities  which  are  necessary  to  use  formula  (2.22) 

(2)  (2) 

are  obtained.  The  second  calculations  for  « , and  a V;  are  made 

n+l  n+l 

with  Simpson's  formula.  An  estimate  of  the  error  is  made  with  expression 
(2.23)  where  mistakes  in  the  numerical  work  or  too  large  a choice  of 
the  increment  produce  a maximum  error  which  is  too  large. 

The  accuracy  of  this  method  is  such  that  larger  increments  can 
be  used  once  enough  starting  values  are  calculated.  Thus,  when  values 
for  »,  a,  and  their  derivatives  are  calculated  up  to  y - 0.3,  enough 
information  is  available  to  use  Milne's  Method  with  increments  of 
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y ■ 0.1.  Table  III  is  constructed  to  show  the  calculations  for  incre- 
ments of  y ■ 0.1. 

Fourth  Step:  Calculation  of  the  outer  boundary  of  the  plate. 

The  boundary  condition  for  w at  the  outer  edge  of  the  plate  i3 
w * n/3  ■ 1.01:72.  In  this  example  the  value  of  u at  y * 1.1  is  greater 
that  n/3. 

An  extrapolation  procedure  is  used  to  obtain  the  outer  boundary. 
The  equation  for  this  extrapolation  is  based  on  Taylor's  expansion 

expressed  in  terras  of  differences. 

. n(n-l)  .2 

“b  “V  ♦-r-t  “ \ 

where  (2.28) 

H ' Vi  • \ • V2  - 2\.i  * v 

The  value  of  "n"  is  2^.ni-3» 

A pictorial  representation  is  given  in  figure  7b.  All  the  quantities 
of  the  above  equation  are  known  except  n.  The  determination  of  n 
gives  sufficient  information  to  obtain  the  outer  boundary  of  the  plate. 
Reorganization  of  formula  (2.28)  gives  this  quadratic  equation  for  n. 

0 • "2  <v?  - 2<w  * * n <-  vs  * K.i  - V - 2 K - V 

(2.29) 

The  quantity  ti^  is  1.0li72.  In  this  example,  the  solution  of  (2.29) 
is  n ■ 2.558.  The  formula  for  the  outer  boundary  of  the  plate  is 

I * yb  * 7k  + n i y.  (2.30) 

In  this  example  - 0.8,  Ay  ■ 0.1,  and  n » 2.558$  therefore,  T is 
T - 1.0558.  Equation  (2.28)  is  reorganized  and  written  in  terns  of  a. 
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% ■ *k  * r <\*2  - 2Vi  * V * I <-  *k*2  * S - 3V  (2-31) 

In  (2.31)  all  the  quantities  are  known  except  a^;  thus  the  value  of 
"a"  at  the  outer  boundary  of  the  plate  can  be  found  with  formula  (2.31). 

Table  III  gives  the  numerical  results  of  thi3  example  while  figure  8 
is  the  graphical  presentation.  With  the  use  of  (2.13),  (2.lli),  and 
(2.19),  equations  (2.1)  are  expressed  as 


ar  m - | ^2  *TT  sin  - n/6) 
«e  “ | ^2  3in  ^ * n/6^  » 


(2.32) 


and  (2.U)  as 


M 


pb^  cos  (to  ♦ n/6). 


(2.33) 


M0 


- I1  - fi.;2.a  >-  pb2  cos  (co  - n/6)  . 

75  l2 


From  the  results  of  the  solution  of  (2.32)  for  a^,  the  deflection 
of  the  mean  surface  of  the  plate,  w(r),  is  calculated  by  a numerical 
integration  of  (2.3U),  the  transformed  equation  I of  Table  I. 


'(r)  ’ f 


roQ  dr  where  w^j  ■ 0. 


(2.3U) 


The  actual  calculations  for  0Q  and  w(r)  are  performed  in  Chapter  6 

w 

where  the  different  approximate  methods  are  compared. 
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TABLE  II 


y 

CO* 

co  - a' 

a 

0 

0 

0 0 

.5 

o.o5 

0.021*8 

0.0006  0.01*28 

0.2*989 

0.10 

o.o5oi 

0.0025  0.0851* 

0.2*957 

o.i5 

0.0766 

0.0056  0.1276 

0.1*901* 

0.20 

0.101*5 

0.0102  0.1692 

0.1*830 

0.25 

0.1353 

0.0161  0.2099 

0.2*735 

0.30 

0.1672 

0.0238  0.21*97 

0.2*620 

TABLE  III 

y 

co’ 

co  - a' 

a 

0.0 

0 

0 0 

0.5 

0.1 

o.o5oi 

0.0025  0.0851* 

0.1*957 

0.2 

0.101*5 

0.0102  0.1692 

0.1*830 

o.3 

0.1672 

0.0238  0.21*97 

0.1*620 

o.i* 

0.22*81* 

0.01*52  0.3257 

0.1*331 

o.5 

0.3678 

0.072*3  0.3952 

0.3970 

0.6 

0.5377 

0.1197  0.1*587 

0.35U5 

0.7 

0.8038 

0.181*1*  0.5103 

0.3055 

0.8 

1.3603 

0.2859  0.51*12 

0.2528 

0.9 

2.1*360 

0.2*733  0.5207 

0.1981* 

0.10 

U.5692 

0.8023  0.3680 

0.1517 

1.0558  1.01*72 


0.1290 
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2.U  Evaluation  of  the  Method  of  Sokolovsky. 

The  method  of  Sokolovsky  pives  a procedure  for  obtaininr  a numerical 
solution  of  equations  II,  III,  and  IV  of  Table  I.  This  solution  is 
obtained  by  the  simultaneous  inte oration  of  the  two  non-linear  first 
order  differential  equations  (2.8)  and  (2.9)«  The  actual  numerical 
integration  of  these  equations  is  greatly  simplified  by  assuming 
expression  (2.11)  for  a stress-strain  relation.  Expression  (2.11) 
is  a parabolic  approximation  for  the  actual  stress-strain  relation- 
ship. The  accuracy  of  the  results  given  in  figure  8 and  9 depends 
upon  the  accuracy  of  this  assumption..  Figure  9 is  a rraph  showing 
the  dependence  of  co  and  *a"  upon  the  radius  of  the  plate,  r,  when 
the  initial  values  of  a ■ 0.25,  a • 0.50  and  a * 1.00  are  used  to 
calculate  the  example  of  section  2.2. 

A check  on  the  accuracy  of  each  step  of  the  numerical  integration 
procedure  is  riven  by  formula  (2.23).  Though  the  error  of  each  step 
is  small,  the  error  of  the  whole  procedure  is  cumulative.  Equations 
(2.32)  and  (2.33)  show  that  the  cumulative  error  in  determining  Y 
alters  the  value  of  the  iioments  and  curvatures  throughout  the  plate. 

When  the  same  size  increment  is  used,  the  behavior  of  the  integration 
process  is  such  that  the  error  is  greater  for  smaller  initial  "a". 

Thus,  this  method  gives  better  results  for  plates  which  are  loaded 
near  to  the  limiting  strain  condition,  or  for  initial  "a"  approxi- 
mately one. 

The  method  of  this  chapter  is  an  inverse  method.  The  variable 
■a",  which  is  a measure  of  the  strain  condition  of  the  plate,  is  riven 
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an  initial  value,  and  the  boundary  condition  (2.19)  determines  the 
parameter  ^ equal  to  \ hpb  /D.  The  parameter  X determines  the 
physical  properties  of  the  plate  and  the  permissible  load.  Any  one  of 
the  factors  in  T can  be  adjusted  to  give  the  correct  boundary  value  of 
T.  But  if  4 plate  of  a given  material  and  specific  dimensions  has  a 
definite  pressure  applied,  T is  determined.  Thus,  a series  of  calcu- 
lations are  necessary  to  adjust  the  initial  value  of  ■a"  so  that  the 
correct  value  of  X occurs  at  the  boundary  of  the  plate. 

The  calculation  of  am  actual  examole,  for  any  given  initial  con- 
ditions, can  be  performed  within  a day.  This  calculation  gives  u> 
and  "a"  as  a function  of  y in  the  form  of  a table  or  a praph.  Thus, 
the  bending  moments,  curvatures,  and  deflections  of  the  plate,  given 
by  (2.32),  (2.33),  and  (2.3U),  cannot  be  expressed  analytically.  A 
complete  description  of  the  physical  behavior  of  the  plate  can  be 
found  to  any  desired  accuracy  when  sufficiently  small  increments  are 
used.  But,  thi3  method  is  too  lengthy  for  quick  estimates  of  the 
moments,  curvatures,  or  displacements  for  plates  with  prescribed 
dimensions.  In  the  chapters  which  follow  several  approximate  methods 
are  described  which  can  provide  these  quick  estimates.  The  last  chapter 
of  thi3  thesis  presents  a comparison  of  the  results  of  these  various 


methods 


CHAPTER  III 


THE  ITERATION  METHOD 


3.1  The  Mathematical  Development  of  the  Method. 

The  paper  by  Ilyushir  (9)  develops  an  approximate  method  for  the 
solution  of  oroblems  with  plastic  stress-strain  relations.  This  method 
leads  to  differential  equations  which  can  be  solved  by  an  iteration 
procedure,  and  is  used  here  to  obtain  an  approximate  solution  for  the 
problem  of  the  bendinr  of  a circular  plate. 

A second  order  ordinary  differential  equation  is  obtained  by 
the  substitution  of  formulae  II  and  IV,  into  III  and  V of  Table  I. 


2 

da  frrp(r)dr  do.  da. 

♦ 'nr ' - i - * B ("  -i  ♦ J t?> 

dr 


{r  XT + 2 ao) 


(3.1) 


where 


o^> 

b -22 

n-l 


o 3X  h2n  f.  /da^2  da  7]  n 

^ |_r  * 3r0©  "3?  + 3 ae  _ 


(3.2) 


The  iteration  method  of  Ilyushin  is  now  applied  to  equation  (3.1). 
Equation  (3.1)  is  reorganized  into 


'1  d<rV 


where 


Qy  ■ p ( rp(r)dr. 


[»  ^ ^ 


b dae 


7 "3 rj 


(3.3) 


(1.5) 
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The  terms  on  the  rirht  iand  side  of  (3*3)  are  separated  into 


7e  ' V 


(3.6) 


and 


E 


(r  -a?  * \ ae> 


b dae) 


* 2 ~Sr 


y 


(3.6) 


The  expression  on  the  left  hand  side  of  (3*3)  is  encountered  in  the 
solution  for  the  elastic  bending  of  a circular  plate.  When  the  stress- 
strain  relation  ic  linear,  the  quantity  B,  given  by  (3.2)  is  zero.  The 
linear  formula  for  a circular  plate  is 

d fl  d / 2 71  Ve 

dr|_r  dr  'r  a0_J  D 

From  (3.U)  and  (3.5),  equation  (3.3)  can  be  rewritten  as 


d 1 1 d <r%>|  - - b <7.  ♦ y • 


3r  I r 3r 


(3.6) 


The  quantity  V , which  is  interpreted  as  a pseudo-shearing  force 
due  to  the  non-linear  stress  strain  relation,  can  be  evaluated  in  terms 
of  an  estimated  Cg  as  a function  of  r.  An  approximate  solution  for 
a_  is  then  found  by  solving  (3.6)  with  the  shear  V and  the  estimated 
p3eudo-shear  V^.  The  approximate  solution  of  this  section  will  not 
be  found  directly  from  (3.3),  but  a convenient  form  suitable  for  an 
iteration  procedure  is  derived  from  (3.3)  by  using  a Green's  Function. 
The  example  of  a simply  supported  plate  is  used  to  demonstrate  this 
method . 

The  boundary  conditions  for  a simply  supported  plate  are: 
a).  At  the  center  of  the  plate 


r - 0,  Mr  - 


Mq,  aQ  is  finite 


(3.7) 


nr"»*  ■'»rr*  ’ 


b).  At  the  outer  edre  of  the  plate 


“Q  1 

r ■ b,  M ■ 0,  r t — * — a.  ■ 0. 
* r * dr  2 9 


(3.8) 


Expressions  (3.7)  and  (3.8)  are  derived  by  substituting  the  given 
values  of  the  moments  into  formulae  III,  Table  I. 

Equation  (3*3)  will  he  expressed  in  an  integral  form  from  which 
Oq  can  be  round  by  integration  for  any  arbitrary  shear  force.  For 
this  integral  representation,  a solution,  satisfying  the  boundary 


conditions  (3.7)  and  (3*8),  is  required  for 


(3.9) 


where 


6 (Jj  ) * 0,  for  0&r<^  £ , ^<r  & b, 

r Pi  d 2 ”1 

and  at  r , 6 (^  ) increases  S d - — (r  aQ)  one  unit.  The  first 
integration  of  (3.9)  yields 


? h <rS>  ■ ci  0ir'? 


(3.10) 


r h (r2“e)  " ci  " 1 ?<rib< 

The  second  integration  of  (3.9)  yields 


C C 

“e  “ "7  * T 
r 


(C1  - 1) 


4 frib- 


(3.11) 


(3.12) 


The  boundary  conditions  (3.7)  and  (3*8)  and  the  condition  that 
is  continuous  at  r ■ J determine  the  three  constants  C^,  C^,  and  Cy 
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These  constants  are 


§ ^ 


C2  ' 


U , C. 
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~T’ 


which  upon  substitution  into  (3.11) 
for  ctg  as 


and  (3.12)  determine  the  solution 


o*  5 

(3.13) 

sb 


The  particular  solution  (3*13)  is  used  to  obtain  the  effect  of 
an  arbitrarily  distributed  shearing  force,  V( J ).  The  magnitude  of 
the  incremental  shearing  force  acting  on  the  area  d$  at  r B f is 


6 (?) 


y(_? ) 

5 


d5  * 


This  shearing  force  produces  an  incremental  change  in  of 
d.9(r)  - 

d oe(r)  - Q-|-)2  * 3 (-f)2]  f<rs.k. 

Due  to  the  linearity  of  equation  (3.9),  the  effects  of  the  incre- 
mental shearing  forces  can  be  integrated  to  form  the  total  increase 
in  Oq  caused  by  the  distributed  shearing  force  V(jf  ).  This  increase 
is 


(3  .ill ) 
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The  integral  form  cf  equation  (3.3)  is  obtained  v\.en  the  expressions 
(3.U)  and  (3.5)  are  sub.°tituted  into  (3.1b)  fcr  the  arbitrary  shearing 
force  V(f). 

V(jT)  - ve  ♦ Vp. 

The  result  of  the  above  procedure,  after  reorganization  ana  several 
integrations  by  parts  to  eliminate  the  derivatives  of  B,  ia 

ae(r)  “ Z I if  (4")2  d ? + I $ if  (4")2  df  * I 5 if  d5 


* 1 f £ (o9 + 4-^r5  M!  ♦ \ f 4 (“9 + 4-ay>  (3.15) 

ar^Bdf* 

The  form  of  this  equation  ensures  the  satisfaction  of  the  boundary 
conditions  for  the  simply  auoported  plate.  Boundary  conditions  (3.7) 
and  (3.8)  are  applied  in  the  formulation  of  the  Green's  Function  which 
determines  for  arbitrarily  distributed  shearing  forces.  For  a 
plate  subjected  to  different  boundary  conditions,  equations  (3.7) 
and  (3.8)  are  replaced  by  new  equations.  The  subsequent  formulae 
are  altered,  but  the  procedure  remains  the  sane. 

Expression  (3.15)  is  the  desired  integral  form  of  (3*3)  for 
the  application  of  the  iteration  method.  The  explanation  of  this 
iteration  process  is  explained  in  the  next  section  by  solving  the 
example  problem. 


i 
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The  integral  form  cf  equation  (3 - 3 ) is  obtained  vhen  the  expressions 

(3.h)  and  (3*5)  are  substituted  into  ( 3 • 3J-* ) for  the  arbitrary  shearing 

force  V(f  ). 

V(  f ) - V + V . 

J e p 

The  result  of  the  above  procedure,  after  reorganization  and  several 
integrations  by  parts  to  eliminate  the  derivatives  of  B,  is 


W) 


5 f r * j f !r  <^>2  df  * I f ir  <*5 


*5  $ ^ (“e  * 4-  ^r>  “S  * I f 4 <«e  * 4-  ay>  “5  o.«) 


The  form  of  this  equation  ensures  the  satisfaction  of  the  boundary 
conditions  for  the  simply  suoported  plate.  Boundary  conditions  (3.7) 
and  (3.8)  are  applied  in  the  formulation  of  the  Green's  Function  which 
determines  Oq  for  arbitrarily  distributed  shearing  forces.  For  a 
plate  subjected  to  different  boundary  conditions,  equations  (3.7) 
and  (3.8)  are  replaced  by  new  equations.  The  subsequent  formulae 
are  altered,  but  the  procedure  remains  the  same. 

Repression  (3.15)  is  the  desired  integral  form  of  (3.3)  for 
the  application  of  the  iteration  method.  The  explanation  of  this 
iteration  process  is  explained  in  the  next  section  by  solving  the 
example  problem. 
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3.2  The  Numerical  Example. 

The  example  oroblem  of  section  2.2  is  used  to  illustrate  the 
iteration  process.  In  this  example  the  load  is  a constant. 
p(r)  - p 


Qr  “ 1/r  J prdr  ■ pr/2 


(3.16) 


A parabolic  stress-strain  relation  i3  assumed  which  determines  the 
quantity  B from  (3.2)  as 


da 


3 V + 3 ae  (r  XT')  * (r  XT') 


da^  2 


B - 0.6  Ah 

With  the  substitution  of  (3-16)  into  equation  (3.15),  a^  is 

[l  - + \[  4 (ae * -r* 


(3.17) 


‘e(r) 


♦ \ r 4 <•« + 4-  i r ^ “F  • 


(3.18) 


o r 


The  iteration  procedure  is  now  applied  to  equation  (3.18). 
a).  The  elastic  solution. 

The  first  step  of  the  iteration  process  is  to  assume  that  A ■ 0. 

This  is  equivalent  to  assuming  that  B is  zero  or  that  the  pseudo- 

shearing  force  7 , equation  (3.5),  is  not  present.  The  aa  calculated 
P H 

from  this  assumption  i3  riven  by  a plate  with  an  elastic  stress-strain 
relation  as 

2 


4e)  ■ & 0 - 3(S>2] . 


which  is  obtained  by  evaluating  (3*18)  when  B - 0. 


(3.19) 
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b).  First  a'rroximation  for  Og. 

The  first  aooroximation  for  Og  is  found  by  integrating  (3.19) 

(e) 

where  (3.18)  is  evaluated  by  substituting  (3.17)  for  B,  and  a'  ' 


from  formula  (3.19)  for  Og.  We  find 

_(1)  _ e . pb2  yU„ 

ae(r)  ae  * CTO  rfi 


V * 

Y2  - X1/2  hpb2/D 


f1  - 0.2327  - 0.287l(|)2  ♦ 0.176U(^)Ll  - 0.0Ul89(£)6. 


(3.20) 


The  dimensionless  parameter  Y is  the  same  parameter  encountered 
in  section  (2.2).  This  step  is  equivalent  to  approximating  the 
pseudo-shearing  force  V , formula  (3.5),  in  terms  of  the  elastic 
effects  on  a0.  The  apparent  shearinp  force  Vp,  due  to  the  plastic 
behavior  of  the  material,  increases  oQ. 

c).  The  second  approximation  for  dg. 

The  value  fc  n. g ' obtained  in  expression  (3.20)  can  be  used  to 
evaluate  B and  oQ  in  (3«17)  and  (3.18)  as  functions  of  r.  Then 
expression  (3.18)  is  integrated  for  the  second  approximation  for 


_ (2)  . (e)  + pb2  Qf 

“q  ae  + UBP  Ly  i 


♦ y f 2 ♦ 


(3.21) 


f2  - 0.0155  - 0.0931  (|)2  ♦ 0.113  (§)U  - 0.0793  (|)6  ♦ 0.0302  (£)8 
- 0.00li5U  (£)10. 

In  this  aDproximation  another  estimation  for  the  apparent  shearing 


force  V is  made  by  using  a rectified  0-.  The  increase  in  the  estimated 
P V 

over  Oq6^  gives  a greater  value  for  a g2^  than  dg^.  The  term 


i 


»r*U  MfeJr 
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0^2  (Y)  denotes  terns  containing  the  twelfth  end  hirher  powers  of  Y. 
These  terns  are  altered  on  further  approximations  for  yet  they 
are  not  needed  for  the  calculation  of  the  twelfth  cower  terns  in  I 


in  the  next  aocroximation. 


d).  The  nth  approxi-ation  for  a 
For  the  n^  approximation  of  o^,  ctg11  ^ i3  used  to  evaluate  the 


G 


intepral  (3.18).  The  form  of  this  solution  is 


>>  - ♦ 


Em  ^ **fk  + °U(n-l)(Y^j  . 


(3.22) 


The  term  T^f  is  the  additional  correction  for  oQ.  This  tern  is  not 
altered  in  further  approximations.  The  terms  of  higher  order  than 
Un  in  the  parameter  Y can  be  neglected  and  do  not  effect  the  next 
approximation  for  the  term  containing  Y^1*^.  The  elimination  of  the 
terms  0 +^(Y)  in  each  approximation  greatly  decreases  the  numerical 

work. 


e).  The  check. 

For  each  approximation  a new  function,  fQ,  is  obtained.  This 

function  is  the  only  term  which  has  the  parameter  Y to  the  lin  cower. 

For  the  computations  to  be  correct,  all  the  functions,  f^,  must 

satisfy  the  boundary  conditions.  After  each  approximation  the  newly 

acquired  function  f should  be  substituted  into  expression  (3.8) 

n 

to  check  the  computations. 

The  functions  f^  and  f^  are  evaluated  in  Table  IV  for  various 

(2 ' 

ratios  of  r/b.  These  values  can  be  used  to  calculate  ' for  various 
loads  and  values  of  the  parameter  Y. 
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TABLE  17 

Evaluation  of  functions  appearing  in 
the  Iteration  Method.  (Formula  3.21) 


r/b 

7 - 3 (r/b)2 

fl  &/»>]] 

f2  [Tr/b^I 

0 

7.00 

0.233 

o.oU55 

0.1 

1 

6.97 

0.230 

0.0iUi6 

m 

6.88 

0.222 

0.0l|20 

0.3 

6.73 

0.208 

0.0380 

o.U 

6.52 

0.191 

0.0332 

0.5 

6.25 

0.171 

0.0282 

0.6 

5.92 

o.i5o 

0.02 3b 

0.7 

5.53 

0.129 

0.0193 

0.8 

5.08 

0.110 

0.0160 

0.9 

U.57 

0.0936 

0.0135 

1.0 

U.oo 

0.0801 

0.0118 

3.3  Evaluation  of  the  Method. 

The  iteration  method  provides  a solution  for  the  plastic 
bending  of  a circular  plate  in  terns  of  the  dimensionless  paraneter  Y 
and  functions  of  r.  In  section  2.2  values  of  Y are  calculated  for 


3U 


various  degrees  of  plastic  yielding.  For  a condition  of  complete 
plastic  yielding  at  the  center  of  the  plate  I ■ 1.368.  Therefore, 
the  ranpe  of  T is  0 — T *-1.386.  The  accuracy  of  this  method  for  a 
given  number  of  aoproximations  depends  upon  the  converrence  of  f 
and  the  value  of  T. 

With  the  use  of  Table  IV,  an  estimate  of  the  maximum  correction 
for  between  a first  and  second  approximation  is  about  7 percent. 

This  estimate  is  obtained  from 

x% 

percent  correction  * 1 100  • (3.23) 

7 - 3(£r  ♦ 

The  error  involved  in  the  two  term  approximation  should  be  less  than 
this  value. 

Each  iteration  increases  the  estimated  value  of  o0,  which  makes 

the  series  of  the  aporoximate  pseudo  shearing  force  a monotonously 

increasing  series.  Physical  considerations  make  it  plausible  that  the 

correct  V is  the  limit  of  this  series  or  that  the  iteration  method 
P 

is  convergent.  Panferov  (5)  treats  the  problen  of  convergence  for 
this  method  mathematically. 

The  exact  solution  for  in  terms  of  the  infinite  series  obtained 
by  the  iteration  procedure  is 

ae  ‘ - 3(E)2  * § ^k]  • 

The  deviation  from  the  nth  approximation  for  the  a ^ given  by  (3.22) 
and  the  exact  is 
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_v2  o*> 

pb  r 1 

^ k-n+1 


(3.25) 


Expression  (3.25)  shows  that  the  sequence  Og  - converges  more 

slowly  the  greater  the  parameter  Y.  This  result  follows  from  the 

consideration  that  the  greater  the  value  of  the  parameter  Y,  which 

indicates  the  extent  of  plastic  straining,  the  greater  is  the  difference 
f 0 } 

a0  - . Thus,  a greater  number  of  approximations  are  needed  to 

correct  the  larger  initial  error  for  large  values  of  the  parameter  Y. 

When  the  desired  accuracy  must  be  preater  than  that  given  by 
a two  term  approximation,  Sokolovsky’s  Method  is  preferable.  The 
numerical  work  for  higher  term  approximations  of  the  iteration  method 
is  excessive.  This  method  has  advantages  over  Sokolovsky's  Method 
because  the  approximations  involve  only  polynomials  in  r which  can 
be  integrated  or  differentiated  to  determine  the  deflection  and  o^. 

The  mathematical  approach  of  the  iteration  method  is  not  as  direct  as 
Sokolovsky's  Method,  but  the  actual  calculations  are  shorter  for  any 
problem  with  a given  load  distribution  and  boundary  conditions  when 
the  two  term  approximation  is  used.  Only  one  set  of  calculations  is 
needed  for  any  outer  radius  of  the  plate  or  magnitude  of  the  para- 
meter determining  the  load  distribution. 


CHAPTffi  IV 


THE  APPLICATION  OF  POTENTIAL  ENHtGY 


U.l  Introduction  to  the  Minimum  Potential  l&iergy  Method. 

When  energy  methods  or  variational  principles  are  used,  the 
variational  principle  should  be  thought  of  as  the  basic  law  govern- 
ing the  behavior  of  the  system.  The  behavior  of  the  system  is  de- 
termined by  a definite  variational  procedure  upon  the  energy  or  equiva- 
lent quantity  of  the  complete  system.  This  method  of  approach  differs 
from  the  procedure  of  the  previous  chapters  where  the  behavior  of 
every  element  is  analyzed.  Both  methods  are  equivalent,  but  the 
assumed  fundamental  concepts  are  different. 

The  variational  principle  used  in  this  c';apter  is  the  principle 
of  minimum  potential  energy.  Ilyushin  (9)  and  Phillips  (10)  give  a 
general  development  of  the  principle  of  minimum  potential  energy. 
Greenberg  (16)  treats  the  subject  of  variational  methods  for  the  theory 
of  plastic  flow  and  the  theory  of  plastic  deformations.  His  report 
gives  a thorough  development  and  includes  the  original  references  for 
the  various  variational  principles. 

The  next  section  directly  approaches  the  problem  of  the  deformation 
of  a circular  plate  using  the  minimum  potential  energy  method.  The 
problem  is  expressed  in  terms  of  the  moments,  curvatures,  and  deflections 
of  the  plate. 

U.2  Formulation  of  the  Minimum  Potential  Energy  Method  for  the  Bending 
of  Circular  Plates. 
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For  the  moment-curvature  relations  IV  of  Table  I,  the  strain 
energy  per  unit  of  a plate  is  expressed  in  terms  of  A as 


H - f ^ . 


(li.l) 


The  potential  energy  for  the  plate  is 


^ ■ 2n  £ (W  - pw)  rdr. 


(U.2) 


For  the  variational  principle  of  this  section  the  following  definitions 
are  necessary* 

1) .  Admissible  displacements  are  displacements  which  agree 
with  the  prescribed  boundary  conditions. 

2) .  Admissible  curvatures  are  curvatures  which  satisfy  equations 
I and  II  of  Table  I throughout  the  plate,  where  the  w in  equation  I 

of  Table  I is  an  admissible  displacement. 

3) .  Admissible  moments  are  moments  derived  from  admissible 
curvatures  by  equations  IV  of  Table  I.  The  load  p is  a given  quantity 
which  is  not  varied.  From  equations  (2.2)  for  A and  (ii.l)  for  W,  we 
find  that ^ can  be  expressed  in  terms  of  only  the  admissible  curvatures 
and  displacements. 

The  principle  of  minimum  potential  energy  states  that  the  admissible 
curvatures  and  displacements  which  render  $ a minimum  yield  admissible 
moments  that  satisfy  the  equilibrium  condition  III  of  Table  I.  Thus, 
these  curvatures  are  a solution  for  the  bending  of  a circular  plate 
which  satisfy  equations  I,  II,  III,  and  IV  of  Table  I.  The  proof  of 
this  principle  follows: 


i 


*?>*•  t 
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A necessary  condition  for  $ to  be  a minimum  is  that  the  first 
variation  of  $ in  zero  for  all  admissible  curvatures  and  displace- 
ments. Hence, 


aA‘ 

35“ 

r 


6a 


♦ B ) “lta 


*9 


rdr 


0 


(U.3) 

From  the  moment-curvature  relations  III  of  Table  I and  formula  (2.2), 
the  moments  are  expressed  in  terras  of  the  curvatures  as 


. P(1  - b)  aA^  . . B(1  - b) 

r 2 55^  * 9 2 35^ 


(U.U) 


The  substitution  of  equations  (U.U)  into  (U.3)  yields 


rb 

) (Mr6ar  ♦ Mg6ag  - p6w)  rdr  - 0 


(U.5) 


The  variations  of  the  displacements  and  curvatures  must  satisfy 
equations  I and  II  of  Table  I: 

bo  - - 

9 r dr  * 


6a  - 
r 


d(r6ag) 

3r 


With  the  substitution  of  the  following  identities 


/b  d(r5aj 

rM  — T_ 

•J  r dr 


dr  - r M 6a_ 
r 9 


- J 


-b  d(rM  ) 


dr 


6aQ  rdr 


(U.6) 


d(rM  ) 


IF 


6aQ  rdr 


(U.7) 
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| 5w  prdr  - 6w  J rpdr  I + ^ ((  rpdr)  [-  i rdr 

o o 'o  o ' o ^ -* 

■ ° * i (f  rpdr-)  E r rdr 


(U.7) 


and  equations  (U.6)  into  (U.5),  we  find 


rb  d(rM  ) r 

) 3r  - M0  M rpdr 

o — o—l 


6a0  rdr  - 0 


(U.8) 


When  expression  (U.8)  is  satisfied  for  all  possible  admi33ible  vari- 
ations of  Oq,  the  bracketed  term  must  be  identically  zero.  This 
bracketed  term  is  the  equilibrium  equation  IV  of  Table  I;  thus,  the 
principle  of  minimum  potential  energy  is  proved. 

U.3  Approximate  Methods. 

Two  approximate  methods  are  available  for  the  solution  of  equations 
(U.2)  and  (U.8).  These  methods  place  restrictions  upon  arbitrary  para- 
meters in  an  assumed  form  of  a solution.  The  explanation  of  these  two 
methods  follows : 

The  first  method,  the  Rayleigh-Ritz  Method,  assumes  an  approxi- 
mate form  of  the  deflection  w with  the  adjustable  parameters  C^, 


(U.9) 


where  each  w^  satisfies  the  boundary  condition  imposed  upon  the  plate. 
An  approximation  for  the  potential  energy  of  the  plate  is  found  by 
using  expression  (U.9)  to  evaluate  the  integral  (U.2). 


I * ( - p(r)w|  rdr 


(U.10) 


Uo 


For  the  calculation  of  W(w),  the  curvatures  are  obtained  with  equations 
I and  II  of  Table  I.  The  parameters  are  varied  to  give  a stationary 

value  of  the  integral  (U.10)  by  setting  its  partial  derivative  with 
respect  to  each  of  the  Darameters  equal  to  zero. 

9I/3C^  - 0 i ■ 1,  n (U.ll) 


Equations  (li.ll)  are  n equations  which  determine  the  n unknown  para- 
meters C^.  The  extremum  of  I is  found  for  the  restricted  set  of 
admissible  displacements  w.  Other  functions  not  contained  in  the 
above  set  may  yield  a lower  potential  energy*  therefore,  this  method 
only  gives  an  approximation  for  the  strain  condition  of  the  plate. 

The  second  method,  Galerkin's  Method,  uses  equation  (U.8)  to 
make  an  approximation  for  oQ.  This  method  assumes  an  approximate 
oQ  with  the  parameters  as 

g ci  “ei(r)-  (‘‘•“l 

Bach  should  satisfy  the  boundary  conditions  of  the  plate.  The 
are  determined  by  evaluating 


C r 

(Qq)  ♦ ] rp(r)dr 

J 


°eirdr  ■ 0;  i 


1, 


n 


(U.13) 

Equations  (U.13)  are  n equations  which  determine  the  parameters  C^. 

Here,  as  in  the  Rayleigh-Ritz  Method,  only  a restricted  set  of  admissible 
curvatures  approximate  <XqJ  therefore,  the  true  strain  condition  is 


not  obtained 


hi 


ELotter  (12)  proves  the  equivalence  of  the  above  two  methods. 

Fox  (13)  (chap.  7)  suggests  what  accuracy  is  obtainable  by  a suitable 
choice  of  assumed  functions.  A general  explanation  of  these  and 
other  approximate  methods  is  given  by  Sokolnikoff  (lh)  (chap.  5). 

The  next  section  demonstrates  the  use  of  Galerkin's  Method, 
h.h  The  Numerical  Example. 

The  numerical  approximation  using  equation  (h.13)  with  Galerkin'3 
Method  is  demonstrated  by  solving  the  example  used  in  sections  2.2 
and  3.2.  This  example  is  a simply  supported  plate  with  a uniform 
load  p and  a parabolic  curvature-moment  relation.  The  determination 
of  the  bracketed  terra  of  formula  (h.13)  as  a function  of  the  assumed 
aQ  is  performed  in  section  3.1.  Therefore,  we  have 


where 


B 


da  da  2~[ 

- 0.6Xh2  3 «Q2  * 3aQ  (r  -~)  * (r  J . 


The  boundary  conditions  are 

a) .  At  the  center  of  the  plate,  r ■ 0, 

M " M_,  and,  a_.  is  finite 
r 6’  ’ 9i 

b) .  At  the  outer  boundary  of  the  plate,  r * b. 


(U.15) 


(U.16) 


da9i  3 

Mr  " r ~3r  + \ a9i  " °* 
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Polynomials  in  r are  selected  as  the  approximating  functions 
The  polynomial  for  a one  tern  anoroxination  is  parabolic.  From 
the  boundary  conditions  xU.16),  oQ1  is  determined  as 


7-3  (r/b)‘ 


(U.17) 


The  form  of  the  elastic  solution  suggests  that  be  written  as 


■ c - 3 <r/»>0 . 


where  C is  the  adjustable  parameter. 


(U.18) 


An  algebraic  equation  for  C is  found  by  evaluating  equations  (li . XU ) 
and  (U.15)  with  expression  (U.18).  This  equation  is 

l\:3  - 0.01925  C ♦ O.OOOU01  - 0,  (U.19) 

where  » X3"^  hpb^/t).  The  parameter  I first  appeared  in  section 
2.2  as  an  indication  of  the  degree  of  plastic  straining. 

Equation  (U.19)  can  be  solved  exactly,  but  the  value  of  C for  a 
given  T is  found  with  less  effort  by  approximate  methods.  The  follow- 
ing method  of  Newton  gives  accurate  results  with  one  approximation. 


Formula  (U.19)  is  expressed  as  a function  of  C. 

f(C)  - rc3  - 0.01925C4-  0.000U01 


(U.20) 


The  (n  ♦ l)th  estimate  of  C is  denoted  in  terms  of  the  n approximation 


plus  a correction  factor  A • 

n 

C n - C ♦ A 
n+1  n n 


(U.21) 


This  correction  factor  is  found  with  the  value  of  the  function  and 

its  derivative  at  C . 

n 


I 


'(CJ 


A - - f(C  ) . 
a n / dC 

The  combination  of  (11.20),  (li.2l),  and  (li.22)  yields 


T^t3  - 0. 


0002005 


^ ' 1.5  ^ - 0. 
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The  values  in  Table  V were  determined  with  formula 


Value  of  the  Multiplvinr  Factor  for 
a One  Tern  Strain-Energy  A oproxi  ration 
TA3L5  V 


T 

C 

0 

.02083 

0.775 

0 .02102 

1.056 

0.02Uj8 

1.368 

0,02305 

For  the  two  term  approximation,  t-ho  nocond  polynomial 
satisfios  the  boundary  condtlona  (b.16),  i* 

&02  * H(r/b)2  - 7(r/b)^. 


The  two  term  approximation  for  ttg  is 


b)J  ♦ 


r 


3(r/br]  ♦ C?  l_ll(r/b)' 


which 


llix 


The  second  term,  a^,  alters  the  va'ue  of  t-e  parameter  found  in 
the  one  term  approximation.  For  the  two  term  esti  ate  of  aQ,  the 
numerical  work  increases  rreatly.  The  form  of  the  algebraic  equations 
for  and  found  by  using  expressions  (h.lii)  and  (U.25)  is 

AUC13  * Wl*  * A13C/C1  * W 4 Vl  4 4 *17  ■ 0 

A21C13  4 Wi*  4 A23C2JC1  4 A2UC23  4 Vl  * A26C2  4 *27  ' 0 

fli.26) 


The  a are  constants  which  contain  the  parameter  T.  For  the  deter- 
mination  of  and  with  a riven  value  of  Y,  two  simultaneous  third 
degree  equations  must  be  solved. 

U.5  Evaluation  of  the  Method. 

A one  term  approximation  by  the  minimum  potential  energy  method 
for  a olastically  bent  plate  adjusts  the  numerical  coefficient  of  the 
elastic  solution.  For  the  elastic  plate  the  value  of  C in  expression 
(h .18)  is  C - l/U8  ■ 0.02083.  Due  to  the  plastic  behavior  of  the 
material,  the  circumferential  curvature  increases  a maximum  of  10.6 
percent.  The  one  term  approximation  gives  a rapid  method  for  directly 
obtaining  a magnitude  correction  factor  for  the  quantities  connected 
with  the  curvature,  but  it  does  not  show  how  the  form  of  the  plate  is 
altered  due  to  the  plastic  behavior  of  the  material. 

The  two  term  approximation  makes  a finer  adjustment  of  the  magni- 
tude factor  and  also  gives  some  information  regarding  variations  of 
the  curvature  with  the  radius.  Besides  the  increased  numerical  work 
necessary  to  obtain  an  expression  of  the  form  (U.26),  two  simultaneous 
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third  degree  algebraic  equations  mu3t  be  solved  for  each  particular 
example.  The  lack  of  a numerical  checking  orocedure  also  presents  a 
strong  objection  to  higher  term  approximations  by  Dotential  energy 
methods.  Further  comparison  of  this  method  with  the  methods  of  the 
previous  chanters  is  riven  in  Chapter  71. 


CHAPTER  V 


THE  APPLICATION  OF  COMPLEMENTARY  POTENTIAL  ENERGY 


5.1  Formulation  of  the  Minimum  Complementary  Potential  Energy  Method 
for  the  Bending  of  Circular  Plates. 

The  material  in  Chapter  IV  introduces  the  variational  procedure. 
The  variational  principle  for  this  chapter  is  the  principle  of  minimum 
complementary  potential  energy.  This  method  is  applicable  when  the 
moment-curvature  relations  state  the  curvatures  in  terms  of  the 
moments.  Thus,  an  inversion  of  equations  IV  of  Table  I is  necessary. 

For  this  inversion  a useful  quantity  M,  similar  in  form  to  the 


stress  invariant  S,  i3  defined  as 

* ■ 5 <•£  - "A  * 0* 


(5.D 


The  substitution  of  equations  IV  of  Table  I into  (5.1)  yields  the 
relation 


(5.2) 


between  A and  M,  the  quantity  A can  be  expressed  in  terms  of  M as 


Tjr  (i  + B’)> 


(5.3) 


where 


»"Cf„  <t>2"  h2°- 

n-1 


(5.U) 


The  Pn's  are  constants  determined  by  substituting  equations  (5*3) 
and  (5.U)  into  IV  and  V of  Table  I and  equating  the  coefficients  of 


like  Dowers  of  M.  The  first  two  coefficients  are 


hi 


h “ | \>  ?2  " 3(|  V?  + 7 V 

Equations  IV  of  Tcble  I,  (5.1),  (5.2),  and  (5.3)  are  3olved  for  the 
curvatures  in  terns  of  the  moments : 

_ 2(1  ♦ B')  3M2 

°r  El  sir  * 

r 


(5.5) 


2(1  + B' 


These  are  the  inverted  moment-curvature  relations. 

For  the  moment-curvature  relatior.3  (5.5),  the  complementary 
potential  energy  per  unit  area  of  the  plate  is  expressed  in  terms 
of  M a3 


wc-( 


(5.6) 


The  complementary  potential  energy  for  the  plate  is 


r b 

) (WC  “ rdr’ 


(5.7) 


For  the  variational  procedure  of  this  section  the  following  defi- 
nitions are  necessary: 

1) .  Admissible  moments  are  moments  which  satisfy  equation  III 
of  Table  I and  the  boundary  conditions  for  the  moments. 

2) .  Admissible  curvatures  are  curvatures  derived  from  admissible 
moments  and  equations  (5.5). 

3) .  Admissible  displacements  are  displacements  calculated  from 
admissible  curvatures  and  equations  I of  Table  I.  Thus,  the  comple- 
mentary potential  energy  (j^can  be  expressed  in  terms  of  only  admissible 


m -i  jx- ic  \ 


U8 


moments  and  the  riven  loads. 

The  princiole  of  minimum  complementary  potential  energy  states 

that  the  admissible  moments  which  render  ® a minimum,  yield  admissible 

c 

curvatures  that  satisfy  equation  II  of  Table  I.  Thus,  the  above 
procedure  gives  the  solution  for  the  bending  of  a circular  plate,  or 
moments  and  curvatures  which  satisfy  equations  I,  II,  III,  IV  of 
Table  I.  The  proof  of  the  principle  of  minimum  complementary  potential 
energy  follows : 

A necessary  condition  for  ^ to  be  a minimum  is  that  the  first 

variation  of  5)  is  zero  for  all  admissible  mo  ients.  Because  the  loads 
Tc 

are  prescribed  on  the  surface  of  the  plate,  the  first  variation  of 


is 


2(1  ♦ B')  3M2  eu 


rdr  - 0. 


(5.8) 


The  substitution  of  equations  (5.5)  into  (5.8)  yields 


) (o6Mr  ♦ agBMg)  rdr  - 0.  (5.9) 

o 

The  variations  of  the  admissible  moments  must  satisfy  equation  III 
of  Table  I: 


d(r5Mr) 


“6M„. 


dr  6 

With  the  substitution  of  the  following  identity 
,b  d(r6Mj 

X6  ”3? 


(5.io) 


/b  d(ron 

1 dr  F “ r^a-6M 


*ev 


tt-f 


■b  d(ra  ) 

-j-H-  6M  rdr 
dr  r 


-f 


d(roe) 

0 - \ - 6M  rdr 

; ar  r 

o 


(5.11) 
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and  equation  (5.10)  into  (5.9)  we  find 

6Mrrdr  - 0.  (5.12) 

When  equation  (5.12)  is  satisfied  for  all  possible  ad  .issible  variations 
of  the  bracketed  tern  nust  be  identically  zero.  Thus,  equation  II 
of  Table  I is  satisfied  and  the  principle  of  minimum  complementary 
energ\  is  moved. 

5.2  Approximate  Methods. 

The  methods  of  section  li.3  can  be  used  to  evaluate  approximately 
the  integrals  (5.7)  and  (5.12).  The  integral  (5.7)  is  bandied  by 
the  Ravleirh-Ritz  Method  with  the  simplification 


f- 


a - 

r 


d(r<IQ) 

3r 


'f  Mc 


rdr  * 0, 


(5.13) 


because  the  pressure  p(r)  is  prescribed  on  the  surface  of  the  plate 

and  does  not  affect  the  variation  of  5 . The  assumed  radial  moment 

c 

M is 
r 


n 


y c m . 

r k in 


(r) 


(5.1b) 


where  the  Mr^(r)  satisfy  the  boundary  conditions  and  the  a~e  the 


adjustable  parameters.  The  circumferential  moment  is 


V 


d (rM  ) 


“3T 


j rPdr, 


(5.15) 


which  assures  the  satisfaction  of  the  equilibrium  conditions  for 
admissible  moments.  The  vanishing  of  the  partial  derivatives  with 
respect  to  the  n parameters  of  expression  (5.13)  evaluated  in 
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terms  of  (5.1b)  and  (5.15)  determines  t'.ie  followin'?  n equations  for 
the  n arbitrary  parameters  C : 


a 

sr 

c o 


} W rdr  = 0.  i ■ 1,  ...,  n 


(5.16) 


Equations  (5.16)  ensure  an  extremium  of  the  complementary  potential 
energy  for  admissible  moment  variations  of  the  restricted  set  (5.1b); 
therefore,  only  an  approximate  solution  is  obtained. 

Equation  (5.12)  is  approximated  by  Galerkin's  Method.  Galerkin's 
Method  gives  n equations  to  determine  the  of  (5.1b)  as 


rf 


« an  - 

r r 


K<\>] 


dr 


M rdr  ■ 0 
ri 


(5.17) 


i “ 1,  . . . , n 


The  curvatures  in  equations  (5.17)  are  calculated  in  terms  of  M with 

r 

the  moment-curvature  relations  (5.5)  and  expressions  (5.1b)  and  (5.15). 
The  result  of  this  calculation  is 

di 2H  dM 


0 -r  T O -a  r 

\ (1  ♦ B')(r  — * 3r  "3JT  + r P + f J rpdr) 

o dr  o 


dM 

* r -3F  (r  IF  * 


y ♦ /*  rpdr)|  M .rdr 
o -I 


(5.18) 


i - 1,  ...,  n. 

These  n equations  are  enough  to  determine  the  n parameters  C^. 

5.3  The  Numerical  Example. 

The  solution  of  the  example  problem  of  sections  2.2,  3*2  and  b.b 
is  now  approximated  by  the  use  of  the  Minimum  Complementary  Potential 
Biergy  Method.  Formula  (5.18)  is  used  for  Galerkin's  approximation 
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method.  This  numerical  example  has  a parabolic  moment-curvature 
curve  when  the  moment  invariant,  M,  is  expressed  in  terms  of  the  curva- 
ture invariant,  A.  The  factor  B'  of  the  inverse  relationship,  becomes 
the  infinite  series 

a/2.  M 2 a/2.  M h 

0.  - 2„U(-  -p— ) ♦ 17.28  ) ♦ ...  (5.19) 


Expression  (5.19)  is  a specialization  of  (5.b)  when  only  X^  is  considered. 
An  approximation  to  expression  (5.19)  of  the  form 


a/2.,,  2 

B»  - 2.1i  (- ■ pi# ) 


(5.20) 


is  made  to  simplify  the  calculations. 

One  justification  of  this  estimate  is  that  the  actual  moment- 
curvature  relationship  of  M to  A is  a parabolic  approximation;  there- 
fore, an  approximation  could  be  made  by  a parabola  of  the  form  (5.20). 
The  factor  X can  be  adjusted  to  give  as  much  accuracy  as  possible  for 
the  range  of  M to  be  considered.  The  adequacy  of  expression  (5.20)  as 
a representation  of  (5.19)  is  discussed  later. 

The  one  term  polynomial  approximation  i3  suggested  by  the  form 
of  the  elastic  solution  of  the  simply  supported  circular  plate  with 
constant  load  p. 


Hr  - Cpb2  Ql  - (r/hf\  . 

Here,  C is  the  adjustable  parameter  and 


(5.21) 


Mpl  - 1 - (r/b)2.  (5.22) 

The  evaluation  of  expression  (5.18)  with  (5.1),  (5.15),  (5.20),  (5.21) 
and  (5.22)  yields  the  cubic  equation  for  C 
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(C3  - 0.7366  C2  + 0.1920  C - 0.0167U) 

♦ 0.5952  C - 0.1302  - 0 , (5.23) 

where  l2  * X3^2  hpb2/D.  The  parameter  I first  appeared  in  section 
2.2  as  an  indication  of  the  degree  of  plastic  straining. 

By  the  approximation  method  of  Newton,  which  is  explained  in 
section  U.U,  the  successive  approximations  for  C are  obtained  from 


n+1 


*k(Cn3  “°*3683  °2  + 0.00837)  ♦ 0.0651 

5 C 2 - 0.7366  C + 0.0960)  ♦ 0.2976 
n 


W. 


(5.21*) 


Fonnula  (5.2U)  i3  u3ed  to  determine  the  values  of  C in  Table  VI. 

Expression  (5.20)  is  now  compared  with  (5.19).  The  quantity 
(\^2M/D)2  is  e valuta  ted  with  (5.1)  and  (5.21)  as 


T 


*1 


:2  [jL  - ii(r/b)2  + 7 (r/b) 
Qr/b)2  - 5(r/b)^J  ♦ (r/b)1*^  . 


a 


(5.25) 


Equation  (5.25)  has  its  maximum  value  at  r * 0. 

(~TT“)2  ^ (5-26) 

This  parameter  indicates  the  difference  between  (5.19)  and  (5.20). 

When  I • 1.368,  the  first  three  terms  of  expression  (5.19)  give  a 
value  of  B'  over  70  percent  greater  than  (5.20).  Thus  for  large  values 
of  T these  two  expressions  are  not  equivalent. 

For  a two  terra  approximation  for  T^,  a similar  procedure  to  that 
of  section  U.ii  can  be  applied.  The  same  numerical  difficulties  arise 
and  the  extra  effort  involved  in  the  calculations  makes  the  two  term 
approximation  impractical. 
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Value  of  the  Multiplying  Factor  for  a One  Term 
Complementary  Strain-Energy  Approximation 

TABLE  VI 


Y 

c 

0 

0.21875 

0.775 

0.2185 

1.056 

0.2178 

1.368 

0.2163 

2.00 

0.209b 

5.U  Evaluation  of  the  Method. 

For  the  calculations  of  section  5.3,  the  moment-curvature  relation 
is  given  by  formulae  (5.2)  and  (5.20).  For  large  values  of  the  oara- 
meter  T,  the  moment-curvature  relation  of  section  (5.3)  is  not  equivalent 
to  the  relation  used  for  the  numerical  calculations  of  the  previous 
chapters.  Thus,  the  numerical  results  of  this  chapter  cannot  be  compared 
with  those  of  the  previous  chapters.  Also,  the  parameter  Y cannot 
be  interpreted  in  terms  of  the  calculations  of  Chapter  II.  With  the  Y and 
C determined  by  equation  (5.23)  and  the  M given  by  (5.25),  H is  a 
monotonously  increasing  function  of  Y;  thus,  Y still  indicates  the 
stress  condition  in  the  plate.  For  the  moment-curvature  relation 
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with  a B'  given  by  (5.20),  the  value  of  Y is  not  restricted,  because 
the  moment-curvature  relation  determined  by  this  B'  does  not  have  a 
horizontal  tangent. 

The  complementary  potential  energy  method  and  the  potential 
energy  method  give  moments  and  curvatures  which  bracket  the  actual 
potential  and  comolementary  potential  energy  of  the  plate.  A discussion 
of  this  bracketing  is  found  in  the  report  by  Greenberg  (l6).  In  this 
work,  the  bracketing  of  the  potential  and  complementary  ootential  energy 
is  not  obtained  because  different  moment-curvature  relations  are  used 
to  calculate  the  potential  and  complementary  potential  energy. 

The  elastic  solution  of  a circular  plate  gives  a where 
C - 7/32  - 0.21875.  This  elastic  coefficient  can  be  compared  with  the 
values  of  C in  Table  VI  which  consider  the  plastic  moment-curvature 
relation.  The  one  term  approximation  of  section  5.3  makes  ar.  adjust- 
ment for  the  magnitude  of  but  does  not  show  the  change  in  the 
distribution  of  the  moments  due  to  plastic  effects.  Higher  term  approxi- 
mations adjust  the  distribution  of  the  moments  but  lead  to  numerical 
difficulties  which  make  them  impractical. 
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CONCLUSION 

6.1  Comparison  and  Interpretation  of  Numerical  Results. 

The  work  of  the  previous  chapters  develops  four  methods  for 
obtaining  apcroximations  of  the  stress  and  strain  condition  of  a 
plastically  bent  circular  plate.  The  numerical  calculations  of  Chapter  V 
for  the  minimum  complementary  potential  energy  method  use  a no  .ent- 
curvature  relation  different  from  that  of  the  previous  chapters. 

This  difference  makes  a comparison  of  the  numerical  results  of  Chapter  V 
with  those  of  the  previous  chapters  meaningless.  Therefore,  only  the 
numerical  results  of  Chapters  II,  III,  and  IV  are  compared. 

The  circumferential  curvature,  O0,  is  u3ed  as  the  measure  of 
comparison  for  the  various  methods.  Figure  10  shows  as  calculated 
by  the  various  methods  for  T M 1,368  and  T - 1.056.  These  cases 
represent  a plate  which  has  reached  a maximum  strain  condition  at  the 
center,  I ■ 1.368,  and  a plate  having  only  half  the  maximum  permissible 
strain  invariant  at  the  center,  Y * 1.056.  The  results  for  Sokolovsky's 
Method  are  calculated  from  the  data  on  Figure  9 and  formula  (2.32). 

The  results  for  the  iteration  method  are  obtained  with  the  use  of 
Table  IV  and  expression  (3*21).  The  curve  for  the  minimum  potential 
energy  method  is  obtained  by  using  Table  V and  formula  (Jj.l8).  The 
iteration  method  is  calculated  by  using  the  two  correction  terms 
f^  and  f^,  while  the  potential  energy  method  uses  a one  term  approximation. 
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For  T ■ 1.056,  the  contribution  of  the  f^  term  in  formula  (2.32) 
of  the  iteration  method  increases  ttg  by  1 percent.  Thus  for  small  Y, 
the  one  term  iteration  approximation  gives  quick,  accurate  re3ult3. 
Actually,  this  approximation  introduces  less  error  for  saalx  Y than 
Sokolovsky’s  Method  when  the  increments  in  the  numerical  integration 
are  of  the  size  used  in  section  2.3.  The  algebraic  work  for  the 
iteration  method  is  also  less  than  that  of  the  minimum  potential  energy 
method. 

The  plate  which  has  reached  the  limit  of  plastic  strain  at  the 
center  offers  the  best  indication  of  the  accuracy  of  the  different 
methods.  The  results  for  this  plate  are  shown  in  figure  10  (Y  ■ 1.368). 
The  values  for  from  Sokolovsky*s  Method  and  the  iteration  method 
compare  favorably.  The  difference  in  these  two  curves  indicates  the 
magnitude  of  the  neglected  terms  in  formula  (3.21).  If  additional 
terms  are  used  in  the  approximation  for  in  the  iteration  method, 

this  value  of  would  approach  the  value  of  oQ  given  by  Sokolovsky's 

Method.  The  values  of  for  the  minimum  -otential  energy  method  lie 
below  those  of  the  other  methods.  When  only  the  first  correction  term, 
f^,  of  formula  (3.21)  is  used,  the  iteration  and  the  minimum  potential 
energy  method  have  close  numerical  agreement. 

The  deflection  of  the  mean  surface  of  the  plate  is  calculated  by 
using  the  results  of  the  iteration  method  in  formula  (2.3U).  The  a g 
from  the  iteration  method  is  U3ed  to  calculate  w(r)  because  it  is 
expressed  in  the  form  of  a polynomial  which  is  easily  integrated  and 
has  a good  numerical  agreement  with  the  given  by  Sokolovsky's  Method. 
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The  result  of  the  substitution  of  (3.21)  into  (2 . 3U ) is 
w(r)-  2.75  - 3.5(r/b)2  ♦ 0.75(r/b)il  ♦ l11  [o.06873 

- 0.1i635(r/b)2  ♦ 0.07l78(r/b)^  - 0.029UO(r/b)^  ♦ 0.00S2U(r/b)®f 

♦ I8  [o.01103  - 0.02275(r/b)2  ♦ 0.02328  (r/b  )**  - O.Ol883(r/b)6 

♦ 0.00991 (r/b)®  - 0.00302 (r/b)10  ♦ C.00038(r/b)^  j (6.1) 

The  values  for  w(r)  in  figure  11  are  calculated  with  formula  (6.1)  and 
give  the  deflection  of  a circular  plate  for  I * 1.368,  I * 1.056,  and 
for  the  case  of  linear  stress-strain  relation  in  which  case  Y ■ 0. 

The  increase  in  the  deflection  at  the  center  of  the  plate  due  to  the 
non-linear  term  in  the  stress-strain  relation  for  Y ■ 1.368  is  about 
13  percent. 

The  radial  bending  moment  M for  Sokolovsky's  Method  is  found  by 

r 

substituting  the  values  for  "a"  and  « of  figure  9 into  equation  (2.33). 
When  the  given  by  formula  (U.18)  is  substituted  into  equation  17  of 
Table  I,  the  for  the  potential  energy  method  is 

Mr  " "h— [l  - (r/b)^  jl  - [39 (r/b )^  - 8U(r/b)2  ♦ iJJ. 

(6.2) 

In  the  example  of  the  simply  supported  plate,  the  boundary  conditions 
on  <Xq  are  determined  by  equation  IV  of  Table  I;  thus,  formula  (6.2) 
gives  an  which  satisfies  the  boundary  condition,  Mr  ■ 0 at  r ■ b. 

In  general,  for  the  minimum  potential  energy  method,  the  admissible 
curvatures  nust  satisfy  the  boundary  conditions,  but  the  derived  bending 
moments  do  not  have  to  satisfy  the  boundary  conditions.  Figure  12 
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gives  a comparison  of  the  radial  bending  moment  calculated  from  equations 
(2.33)  and  (6.2)  for  the  fully  strained  plate,  where  Y ■ 1.368. 

6.2  Remarks  on  the  Theory  of  Plastic  Flow. 

In  section  l.h  the  stress-strain  relations  of  the  secant  modulus 
theory  (the  theory  of  plastic  deformations)  are  formulated.  The 
combination  of  (1.13)  and  (l.lli)  gives 

Si  * 20  H ' f(E"3  V 1 " 1'2’3*  (6.3) 


Inverse  relations  relating  the  total  strains  to  the  stresses  are  found 
by  a orocedure  similar  to  that  of  section  5.1.  These  inverse  relations 
are  of  the  form 


e 


i 


s±,  i - 1,2,3. 


(6.U) 


Repressions  (6.3)  and  (6.U)  are  valid  for  the  case  of  loading  only. 
They  relate  the  final  stress  condition  to  the  final  strain  condition 
of  the  element  independently  of  the  path  of  loading. 

Recent  experimental  investigations  by  Phillips  (20)  and  Morrison 
and  Shepard  (21)  indicate  that  the  plastic  behavior  of  a metal  is 
described  by  the  laws  of  plastic  flow.  Hill  (7)  (chap.  2)  and  Prager 
and  Hodge  (18)  (chap.  1)  give  a formulation  of  the  theory  of  plastic 
flow.  In  what  follows  we  shall  consider  the  case  of  loading  only. 

The  theory  of  plastic  flow  relates  the  increments  of  the  strain 
deviation  to  l).  the  increments  of  the  stress  deviation  and  2).  the 
stress  invariant  by 

de±  - ♦ SjhOJdsj  , i - 1,2,3.  (6.5) 
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Thus  the  total  strain  deviations  are  given  by  the  line  integrals 


1 


<>.) 


s^CsJdS 


, i - 1,2,3, 


(6.6) 


which  depend  upon  the  manner  of  stressing  the  element.  The  appear- 
ance of  the  line  integrals  in  (6.6)  makes  the  inversion  of  the  stress- 
strain  relations  impossible.  Hill  (7)  (chap.  2)  shows  how  to  exore3S 
the  increments  of  the  stress  deviations  in  terms  of  the  stresses  and 
the  increments  of  the  strain  deviations  as 


, i “ 1,2,3.  (6.7) 

Equation  (6.7)  does  not  determine  the  increments  of  the  stress  deviations 
in  terms  of  the  strains  only,  but  also  assumes  knowledge  of  the  stress 
condition  of  the  element. 

Because  the  stress  condition  of  the  body  must  be  known  to  determine 
increments  in  the  stresses  for  given  increments  in  the  strains,  a 
solution  for  the  bending  of  a plate  which  satisfies  assumption  c). 
section  1.3,  must  be  a step  by  step  procedure.  For  this  procedure, 
increments  in  the  curvatures  determine  increments  in  the  strain 
deviations  which  satisfy  assumption  c).  These  strain  deviations 
together  with  the  knowledge  of  the  stress  condition  from  the  previous 
step,  can  be  used  to  calculate  the  increments  of  the  stresses  from 
equation  (6.7).  These  stress  increments  are  used  to  determine  incre- 
ments of  the  bending  moments  which  must  satisfy  the  equilibrium  equation 
III  of  Table  I.  The  numerical  work  of  such  a procedure  is  enormous . 


ds. 


2G 


/^(S)  ^ Sjdej^ 

dei  V 
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An  estimate  is  now  made  of  the  change  in  the  numerical  results 
of  the  previous  chapters  if  the  laws  of  plastic  flow  were  used  instead 
of  the  laws  of  plastic  deformations.  For  this  estimate,  figure  13 
is  constructed  by  the  procedure  similar  to  that  of  section  2.1  for 
the  construction  of  figure  5b.  Then,  figure  13  is  used  as  an  aid  to 
differentiate  between  the  laws  of  plastic  flow  and  plastic  deformations. 

For  the  construction  of  figure  13,  the  strain  deviations  are 
divided  into 

e±  - i - 1,2,3,  (6.8) 

where  e^®  is  the  elastic  strain  deviation  given  by 

e^  - 8i/2G  i - 1,2,3,  (6.9) 

for  both  (6.U)  and  (6.5).  The  plastic  strain  deviation  e^  is  given 
by 

n a<s(s) 

ei  * "TB — 1 “ V»3,  (6.10) 

for  the  theory  of  plastic  deformations,  and  by 
s h(S)  dS 

d0iP  " — 1 " V.3,  (6.11) 

0 

for  the  theory  of  plastic  flow.  Because  e^  is  the  3ame  for  both 
theories,  only  e^*5  is  compared  to  determine  the  difference  between  the 
theories  of  plastic  deformations  and  plastic  flow. 

In  figure  13a,  axes  (VT  e^/2),  (-/3  «2p/2),  and  (VT  e-jP/2 ) lie 
in  the  principal  stress  axes  <T^,  and  <Ty  The  angle  u is  determined 
by  a vector  where"? has  the  components  s^/vT,  and  s^A/T. 

The  vector  ? lies  in  the  octahedral  plane  and  is  equal  in  magnitude 
to  the  stress  invariant  S. 


\ 


t 

f 

[ 

I 

The  surface  shown  in  figure  13b  is  a surface  of  revolution  about 
the  S axis  which  determine®  the  strain  hardening  characteristics  of 
the  material.  The  S axis  is  perpendicular  to  the  octahedral  plane 
described  above.  For  the  theory  of  plastic  deformations,  the  value 
of  Kp 

_ 3 1/2 

eP  - (1/2  e/6^)  (6.12) 

is  the  radius  of  the  circle  formed  by  the  intersection  of  a plane 
parallel  to  the  octahedral  plane  and  a distance  S from  the  octahedral 
plane.  The  substitution  of  equations  (6.10)  into  (6.12)  yields 

gP  . , (6.13) 

where  equation  (6.13)  determines  the  curve  which  generates  the  strain 
hardening  surface.  The  orthogonal  projections  of  a vector  Ep,  at  the 
angle  to  in  the  octahedral  plane,  on  the  axes,  ( -/J  e^p/2),  ( VT  e„p/2), 
and  ( V3  e^p/2 ) are  the  plastic  strain  deviations  given  by  equations 
(6.10). 

For  the  theory  of  plastic  flow,  expressions  (6.1l)  are  used  to 
formulate  the  quantity 

dS.  (6.11*) 

This  quantity  is  the  increment  in  vhe  radii  of  the  circles  formed  by 
the  intersections  of  the  planes  parallel  to  the  octahedral  plane  at 
the  distances  S and  S ♦ dS  with  the  strain  hardening  surface.  Equation 
(6.1U)  determines  the  shape  of  the  curve  which  generates  the  strain 
hardening  surface.  The  orthogonal  projections  of  the  incremental  vector 


f 


i 
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/ J 

/y'  1/2  } ^ de^de^,  at  the  anrle  w in  the  octahedral  plane,  on  the 

axes  (V3  e^/2),  ( */3  e2p/2),  and  (VJ  e^p/2)  are  the  increments 
of  plastic  strain  deviations  siren  by  equations  (6.11). 

As  a simplification,  the  principal  axes  are  assumed  to  reiain 
fixed  throughout  loading.  This  assumption  is  true  for  the  case  of 
circular  plates  symmetrically  loaded.  Also,  the  strain  hardening 
surface  is  assumed  the  sa  is  for  the  theory  of  plastic  flow  and  for  the 
theory  of  plastic  deformations.  With  these  assumptions,  and  if  u ■ const, 
during  loading,  (increasing  S ) , the  same  total  plastic  strain  deviations 
are  predicted,  by  both  of  the  above  theories.  But  when  « changes  during 
loading,  the  value  of  the  total  plastic  strain  deviations  given  by  the 
theory  of  flow  will  depend  upon  how  u>  changes,  while  the  total  plastic 
strain  deviations  given  by  the  theory  of  plastic  deformations  depend 
upon  the  final  value  of  u.  Thus,  for  any  case  of  loading  of  an  element 
where  the  anrle  u>  changes  during  loading,  the  values  of  the  total  plastic 
strain  deviations  predicted  by  the  theory  of  plastic  deformations  differ 
from  those  predicted  by  the  theory  of  plastic  flow.  However,  when  the 
total  change  of  w is  slight,  this  difference  is  small  enough  to  be 
neglected. 

The  variable  u introduced  in  section  2.1  was  determined  by  the 
direction  of  a vector  IT  in  the  octahedral  plane.  Because  the  components 
of  are  proportional  to  the  components  of  E*,  the  vector  A has  the  same 
direction  as  the  vector  E.  This  vector  E is  parallel  to  the  vector  S 
for  the  laws  of  plastic  deformations;  thus,  the  variable  w used  in 
section  2.1  is  the  same  as  the  above  «.  An  inspection  of  figure  9 


.avs 
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shows  that  the  chanre  in  cc  is  only  a few  degrees  at  any  given  value 
of  r/b  during  the  loading  of  the  plate;  hence,  the  moments  predicted 
by  the  theory  of  plastic  deformations  should  be  very  close  to  those 
predicted  by  the  theory  of  plastic  flow,  but  not  identical. 

6.3  The  Application  of  the  Methods  of  Chapters  III,  IV,  and  V to 
the  Bending  of  Rectangular  Plates. 

The  small  deflection  theory  of  rectangular  plates  can  be  treated 
with  the  methods  of  Chapters  III,  IV,  and  V.  A brief  outline  of  the 
application  of  these  methods  is  now  given.  Swida  (19)  develops  the 
kinematic  and  equilibrium  equations.  Ilyushin  (9)  applies  the  method 
of  minimum  potential  energy  and  the  iteration  method  to  the  problem 


* 


6U 


In  the  rectangular  plate  theory  e , e , and  <r  are  assumed  to 

xz  yi  z 

be  zero.  For  the  following  formulation,  the  material  is  assumed 
incompressible,  that  is 


e ♦ e ♦ e - 0. 
xx  yy  zz 

From  the  abov9  assumptions,  the  stresses  are  written  in  terms  of 
the  strains  as 


(6.17) 


<rx  - liOp.  - F(E)J  (e^  ♦ *^2) 
ay  - UQ  [l  - F(iJ  («yy  ♦ e^/2) 

20  (I  - F(b7|  e . 


y 

r 

*y 


(6.18) 


The  second  strain  invariant  B is 

E2  - i J ] e e - e 2 ♦ e 2 ♦ e e ♦ e 2, 
2 “ ij  ij  xx  yy  xx  yy  xy 


(6.19) 


and  as  in  Chapter  II,  E can  be  written  in  terms  of  a quantity  A as 

E2  - z2A2  (6.20) 

where 

A2  - w 2 ♦ w 2 ♦ w w ♦ w ‘ 2.  (6.21) 

xx  yy  xx  yy  xy 

The  function  F(E),  which  expresses  the  plastic  properties  of  the  material. 


is  given  by  the  polynomial 


F(B)  - 


E211  - 


v — ■«  2n.2n 

> X z A . 
* v n 
n*l 


n-1 


(6.22) 


The  definitions  of  the  bending  moments 
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l V**  * V| 


O'  zdz 


, M - $ 

9 TV  ' 


xy 


T zdz 

xy 


(6.23) 


I 


65 


are  combined  with  equations  (6.16),  (6.18),  and  (6.22)  to  yield 
- - 0(1  - B)  («„  ♦ v^), 

Hy  - - D(1  - B)  (.^  ♦ *J2), 

«Iy  * “0(1  “ B)  .^2, 


whers 


D - 8G  h3/3 


<y  31 


'k 


(hk)dn. 


The  equ librium  equation  for  a rectangular  Dlate  is 
02M  32H  d2M 

— jp  + 2 + — 2^  + P<x»y)  " °t  (6.26; 

dx  3 dy 

where  p(x,y)  is  the  distributed  load  per  unit  area  which  acts  normal 
to  the  surface  of  the  plate.  Equations  (6.16b),  (6.2U),  and  (6.26) 
are  the  equations  which  are  used  to  find  the  solution  for  the  bending 
of  a rectangular  plate. 

A formula  which  can  be  used  for  the  iteration  method  is  obtained 
by  the  substitution  of  equations  (6.16)  and  (6.2ii)  into  (6.26).  These 
substitutions  yield 


vllv  " p (p(»?y)  + 


where 


!i  Jx 

* + 244 


dx2dy^  dj* 
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and 


p - D B 
P 


li  , * 32B 

♦ t,w  ♦ w ,„)  — *•  ♦ w 
XX  yy/2  ^2  xy 
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♦ (w  ♦ v 

yy 


v d2B 

W2> 


(6.28) 


pp  is  a pseudo-loading  due  to  the  non-linearity  of  the  stress-strain 
relations.  If  the  iteration  method  of  Chapter  III  were  followed,  a 
Green's  Function  for  the  elastic  plate  would  be  developed  so  the 
effect  of  any  arbitrary  loading  on  the  deflection  w could  be  found  by 
evaluating  an  integral.  In  order  to  find  this  Green's  Function,  the 
elastic  plate  problem  of  a concentrated  load  acting  at  any  arbitrary 
point  on  the  plate  with  the  given  boundary  conditions  must  be  solved. 

Such  solutions  are  given  in  terms  of  infinite  series  which  make  the 
evaluation  of  very  lengthy.  A direct  approach,  which  involves  less 
numerical  effort,  is  to  first  assume  B * 0 and  directly  solve  the  elastic 
problem 


k,  - 


(6.29) 


for  the  given  boundary  conditions.  The  w , found  above,  is  then  used 
to  evaluate  p^  with  (6.28).  This  Pp  is  used  to  obtain  a first  approxi- 
mation for  w with 


^7^  ■ p (p(x,y)  ♦ pp).  (6.30) 

In  general,  analytic  solutions  of  (6.29)  and  (6.30)  are  difficult  to 
find,  but  relaxation  or  methods  of  finite  differences  can  be  used  to 
get  numerical  results.  The  w found  above  can  be  used  to  evaluate  a 
new  p for  a second  approximation  of  w.  Judging  from  the  numerical 
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results  found  in  Chapter  III,  only  two  iterations  should  give  sufficient 
accuracy. 

The  minimum  potential  energy  method  of  Chapter  17  can  be  used  to 
obtain  an  approximate  solution  for  the  bending  of  a rectangular  plate. 
From  the  combination  of  (6.21)  end  (6.2li),  we  find  the  moment-curvature 
relations  in  terms  of  A as 


M - 
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D(1  - B) 
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a2A2 
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xy 


(6.31) 


For  the  moment-curvature  relations  (6.3l)>  the  strain  energy  per  unit 
area  of  the  plate  is 


o 

The  potential  energy  for  the  plate  is 


(6.32) 


$ - Igl  (W  - pw)dxdy.  (6.33) 

The  first  variation  of  $ must  be  zero  for  minimum  potential  energy.  When 
equations  (6.31)  are  substituted  into  6|j  * 0,  the  result  is 

- 4^  Mx^”6wxx^  + 2Mxy^wxy^  + “ P6"  dxdy  ” 0 (6.3U) 

Expression  (6.3k)  can  be  put  in  an  alternate  form  which  shows  that 
the  minimization  of  the  potential  energy  for  admissible  curvatures 
satisfying  (6.l6b)  and  deflections  satisfying  the  given  boundary 
conditions  yield  moments  which  fulfill  the  equilibrium  conditions. 
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To  do  this,  integrations  by  parts  of  the  type 


■ JLx** ds  * y ^ 6 * ^ 


C aM  /-  r 

0 O ^ 6-da  - j j 


a2M 


Boundary 


S ax 


x . d2M 

6wdxdy  - 0 - j ( — Bwdxdy 
S dx 


(6.35) 


are  performed.  Then  equation  (6.3U)  becomes 


' d2M 


32M 


a2M 


(6.36) 


's'  \ax‘ 

The  bracketed  term  of  (6.36)  is  the  equilibrium  condition  (6.26). 

An  approximate  solution  is  obtained  by  assuming  a form  of  the 
deflection 
n 


* ■ ZZ  c<vA*>y) 

i-l  1 1 


(6.37) 


where  the  are  adjustable  parameters  and  the  w^  satisfy  the  given 
boundary  conditions  of  w.  The  can  be  evaluated  by  the  Rayle igh-Ritz 
method  where  $ i3  evaluated  in  terms  of  w and  the  n equations 

0 i-l,  ...,  n (6.38) 


ai 

*7 


determine  the  C^.  The  equivalent  method  of  Galerkin  may  be  used  with 
equation  (6.36). 

The  minimum  complementary  potential  energy  -.ethod  cf  Chapter  7 
can  be  applied  to  a rectangular  plate  for  approximate  solutions. 
Expressions  (6.2li)  determine  a quantity  M as 


M2  » i (M  2 ♦ M 2 - M M ♦ 3M  2) 
3 x y x y ^ xy 


■ [0.11-31 
_ 2 _ 


(6.39) 
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As  in  section  (5.1),  A is  expressed  in  terms  of  M as 
A “ | (1  ♦ B')M 


where 


B«  - 


B (S£)2n  M2" 


(6.1*0) 


(6.1*1) 


The  B have  the  same  values  as  in  section  5.1.  Inverse  moment-curvature 
rn 

relations  are  obtained  from  (6.2l*),  (6.39),  and  (6.1*0)  as 
1*  tn \ _ 2(1  ♦ B')  3M2 

rnry  <MX  - y2)  — ^ — L m-  > 


■ Wxx  " TO" 

Jk. 

...  U txt  u _ 2(1  ♦ B»)  dM2 
" Wyy  3Kl“- (My  V2)  3J T ’ 


(6.1*2) 
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xy 


xy  xy 

For  the  moment-curvature  relations  (6.1*2),  the  complementary  strain 
energy  per  unit  area  of  the  plate  is 


. f 


(6.1*3) 


and  the  complementary  potential  energy  of  the  plate  is 

j>c  ■ y ^ (Wc  - pw)dxdy.  (6.1*1*) 

S 

The  first  variation  of  $c  must  be  zero  for  minimum  complementary 
potential  energy. 

For  the  computational  work,  both  M and  M can  be  assumed  as 

x y 


\ g CiHxl  <*#)  > V ^ Vji  {l’7h 


(6,1.5) 
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where  and  are  parameters,  and  and  satisfy  the  given 
boundary  conditions  on  the  moments.  The  twistinp  moment  is  calculated 
from  (6.26)  as 


V ' ? Oftfr*5  ***  * 5? Ci  ^ dy 


(6.U6) 


♦ F(x)  ♦ G(y), 


where  F(x)  and  G(y)  are  chosen  so  that  they  give  the  correct  boundary 

conditions  for  M . When  the  load  is  prescribed,  gjy  is  zero;  thus, 

the  term  vp  can  be  disregarded  in  the  calculation  of  $ because  it 

^ c 

has  no  effect  on  the  first  variation  of  $ . i>  is  calculated  with 

c c 

equations  (6.39),  (6.1*3),  (6.1*1*),  (6.1*5)  and  (6.1*6)}  then  the  and 
can  be  obtained  from  the  2n  equations 


and 


(6.1*7) 


This  section  has  developed  methods  for  the  approximate  solution 
of  the  small  deflection  theory  of  rectanpular  plates  which  parallel 
those  for  the  small  deflection  theory  of  circular  plates.  Similar 
solutions  for  the  large  deflection  theory  of  plates,  where  the  stretch- 
ing of  the  neutral  surface  is  considered  are  not  readily  obtainable. 
Even  the  elastic  solution  of  such  problems  gives  non-linear  equations 
which  are  difficult  to  solve,  and  as  we  have  seen  the  elastic  solution 
is  a first  step  in  the  iteration  process.  In  addition,  due  to  the 
terms  from  the  stretching  of  the  neutral  surface,  the  strain  invariant 
E cannot  be  expressed  as  a function  of  the  curvatures  only.  Therefore, 
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the  moments  and  curvatures  cannot  be  the  only  variables  in  the 
formulization  of  the  differential  equations.  Because  of  these 
mathematical  difficulties,  the  condition  of  a plate  is  usually 
approximated  by  either  the  small  deflection  theory  of  bending, 
or  a large  deflection  theory  where  only  membrane  stresses  and 
strains  are  considered. 
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Variation  of  *a"andW  with 
Plastic  Banding  of  a Simply 
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